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ABSTRACT 

The basic mechanisms responsible for producing winds from cool, late-type stars are still largely unknown. 
We take inspiration from recent progress in understanding solar wind acceleration to develop a physically 
motivated model of the time-steady mass loss rates of cool main-sequence stars and evolved giants. This 
model follows the energy flux of magnetohydrodynamic turbulence from a subsurface convection zone to its 
eventual dissipation and escape through open magnetic flux tubes. We show how Alfven waves and turbulence 
can produce winds in either a hot corona or a cool extended chromosphere, and we specify the conditions 
that determine whether or not coronal heating occurs. These models do not utilize arbitrary normalization 
factors, but instead predict the mass loss rate directly from a star's fundamental properties. We take account of 
stellar magnetic activity by extending standard age-activity-rotation indicators to include the evolution of the 
filling factor of strong photospheric magnetic fields. We compared the predicted mass loss rates with observed 
values for 47 stars and found significantly better agreement than was obtained from the popular scaling laws 
of Reimers, Schroder, and Cuntz. The algorithm used to compute cool-star mass loss rates is provided as 
a self-contained and efficient computer code. We anticipate that the results from this kind of model can be 
incorporated straightforwardly into stellar evolution calculations and population synthesis techniques. 
Subject headings: stars: coronae — stars: late-type — stars: magnetic field — stars: mass loss — stars: winds, 
outflows — turbulence 



1. INTRODUCTION 

All stars are believed to possess expanding outer atmo- 
spheres known as stellar winds. Continual mass loss has 
a significant impact on the evolution of the stars them- 
selves, on surrounding planetary systems, and on the evolu- 
tion of gas and dust in galaxies (see review s bv lDupreelll986t 
lLamers & Cassinellill999l:lPuls et alj2008l) . For example, the 
Sun's own mass loss was probably an important factor in the 
early erosion of atmospheres from the inne r planets of our so- 
lar system (e.g., Wood 2006; Glide] 2007). On the opposite 
end of the distance scale, a better understanding of the winds 
from supergiant stars is leading to new ways of using them as 
"standard candle s" to measure the distances to other galaxies 
(Kudritzki 2010). By studying the physical mechanisms that 
drive stellar winds, as well as their interaction with processes 
occurring inside the stars (convection, pulsation, rotation, and 
magnetic fields), we are able to make better quantitative pre- 
dictions about a wide range of astrophysical environments. 

Over the last half-century, there has been a great deal 
of research into possible mechanisms for driving stellar 
winds on the "cool side" of the Hertzsprung-Russell di- 
agr am; i.e., effective tem p eratures less than about 800 
K dHolzer & Axfordl 119701: lHartmann & MacGregol [198 
Hear nlll98a lLafon & Berruveril 1991b lMulladll99ft IWillso: 
2000; Hol zwarth & Jardinal2007l) . Despite this work, there 
is still no agreement about the fundamental mechanisms re- 
sponsible for producing these winds. Many studies of stellar 
evolution use approximate prescriptions for mass loss that do 
not depen d on a true physical model of how the outflow is 
produced dReimersll 1 9751 ; lLeithererl2010 ) . Observational val- 
idation of models is made difficult because mass loss rates 
similar to that of the solar wind (M ~ 10~ 14 Mq yr" 1 ) tend to 
be too low to be detectable in most observational diagnostics. 

Fortunately, there has been a great deal of recent 



progress toward identifying and characterizing the pro- 
cesses that produce our own Sun's wind. Self-consistent 
models of turbulence-driven coronal heating and solar 
wind acceleration have begun to succeed in reproducing a 
wide range of observati ons without the need for ad hoc 
;e parameter s (e.g. , ISuzukil l2006t ICranmer et"aT] 120071; 
Rappa zzcTet all 120081 IVerdini et al.l l2010t iBingert & PeterJ 
201 U Ivan B allegooii en et alj|201 ll) . This progress on the so- 
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lar front provides a fruitful opportunity to better understand 
the fundamental physics of coronal heating and wind acceler- 
ation in other kinds of stars. 

The goal of this paper is to construct self-consistent physi- 
cal models of cool-star wind acceleration. These models pre- 
dict stellar mass loss rates without the need for observation- 
ally constrained normalization parameters, artificial heating 
functions, or imposed damping lengths for waves. We aim to 
describe time-steady mass outflows from main-sequence stars 
with solar-type coronae and from giants with cooler outer at- 
mospheres. In princi ple, then, these models cross the well- 
known dividing line (iLinskv & Haischlll979l) between stars 
with and without X-ray emission. However, there are sev- 
eral types of late-type stellar winds that our models do not at- 
tempt to explain: (1) Highly evolved supergiants and asymp- 
totic giant branch (AGB) stars presu mably have winds driven 
by radiation p ressure on dust grains ( lLafon &Berruveill 99 U 
lHofneril20Tll) and/or strong radial pulsations (IWillsonll2000l) . 
(2) T Tauri stars have polar outflows that may be energized by 
magnetosp heric streams o f infalling gas from their accretion 
disks (e.g., Cranmer 2008). (3) Blue horizontal branch stars 
may have line-d riven stellar win ds similar to those of O, B, 
and A type stars dVink & Cassisill2002l) . 

The remainder of this paper is organized as follows. In 
Section [2] we outline the relevant properties of Alfven waves 
and magnetohydrodynamic (MHD) turbulence that we expect 
to find in cool-star atmospheres. Section [3] presents deriva- 
tions of two complementary models of mass loss for stars with 
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and without hot coronae, and also describes how we estimate 
the total mass loss due to both gas pressure and wave pres- 
sure gradients. Section [4] summarizes how we determine the 
level of magnetic activity in a star based on its rotation rate 
and other fundamental parameters. We then give the result- 
ing predictions for mass loss rates of cool stars in Section [5] 
and compare the predictions with existing observational con- 
straints. Finally, Section [6] concludes this paper with a brief 
summary of the major results, a discussion of some of the 
broader implications of this work, and suggestions for future 
improvements. 



2. ALFVEN WAVES IN STELLAR ATMOSPHERES 



been 



For several decades, MHD fluctuations have 
studied as likely sources of energy and momen- 
tum for accel e rating winds from cool stars (see, e.g., 
Hollwed 119781: lHar tmann & Mac Gregorl Il980t IPeC ampli 
198U Iwang & Sheelevl 119911: lAirapetian et al.1 l2000t 
Falceta-Goncalves et al.l l2006t ISuzukil 12007b . Specifically, 
the dissipation of MHD turbulence as a potential sourc e of 
heating for the solar w ind goes back to IColemanl d 1968b and 
Jokipii & Davis d!969l) . Despite the fact that other sources of 
heating and acceleration may exist, we choose to explore how 
much can be explained by restricting ourselves to just this 
one set of processes. The ideas outlined here will be applied 
to both the "hot" and "cold" models for mass loss described 
in Section[3] 

2.1. Setting the Photospheric Properties 

We begin with five fundamental parameters that are as- 
sumed to determine (nearly) all of the other relevant proper- 
ties of a star: mass M*, radius R*, bolometric luminosity L„, 
rotation period P m , and metallicity. We also assume that the 
star's iron abundance, expressed logarithmically with respect 
to hydrogen, is a good enough proxy for the abundances of 
other elements heavier than helium; i.e., [Fe/H] « log(Z/Z Q ). 
For spherical stars, the effective temperature T e g and surface 
gravity g are calculated straightforwardly from 



aT eff = 



AttRI 



GM* 
Rl 



(1) 



where a is the Stefan-Boltzmann constant and G is the New- 
tonian gravitation constant. 

We need to know the mass density in the stellar photo- 
sphere p„ in order to specify the properties of MHD waves 
at that height. The density was computed from the criterion 
that the Rosseland mean optical depth should have a value 
of 2/3 in the photosphere. We used the AESOPUS opacity 
database, 1 which is a tabulation of the Rosseland mean opac- 
ity kr as a function of temp erature, density, and metallicity 
(see lMarigo & Aringerl2009r) . An approximate expression for 
the Rosseland optical depth tr in the photosphere, 



tr = K^pM* = 2/3 



(2) 



was solved for p*, where H* is the photospheric value of the 
density scale height (see below). We used straightforward lin- 
ear interpolation to locate the relevant solutions for p* as a 
function of T e ff, g, and [Fe/H] . 

Figure [Tf a) shows how the photospheric density varies as a 
function of T e ff and log g under the assumption of solar metal- 
licity ([Fe/H] = 0). For context we also show the location 

' |http : //stev . oapd. inaf . it/cgi-bin/aesopus| 
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FIG . 1 . — Derived photospheric parameters shown as a function of both T c ff 
and log g. Contour labels denote (a) base- 10 logarithm of mass density in 
g cm" 3 , (b) magnetic field strength in G, and (c) Alfven wave amplitude 
Vx» in km s~' . Also shown is a post-main-sequence evolutionary track for a 
1 Mq star (red curves) and the location of the ZAMS (white and gray curves). 

of the zero-age mai n sequence (ZAMS) from the models of 
iGirardi et alJ (l2000h . as well as a post-main-sequence evolu- 
tionary track for a 1 star from the BaSTI 2 model database 
dPietrinferni et alj|2004 . 

We used the 2005 release of OPAL plas ma equations of 
state 3 (see also iRogers & Nayfonovl [2002) to estimate the 
mean atomic weight p in a partially ionized photosphere. For 
the range of parameters appropriate for cool stars, we found 
that p is primarily sensitive to r e ff , and not to gravity or metal- 
licity, so we produced a single parameter fit, 



7 1 u /3500-r eff 

u w - + - tanh 

P 4 2 V 600 



(3) 



where T e s is expressed in K. Other quantities that will be 
needed later include the photospheric density scale height, 
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which is given by 



H, = ^ (4) 

where k% is Boltzmann's constant and mu is the mass of a 
hydrogen atom. We also need to compute the equipartition 
magnetic field strength, 



B eq = y/SnP, = 



l 8npi,k B T ef f 



(5) 



where P* is the photospheric gas pressure. Because T e ff and 
\i do not vary over many orders of magnitude, it is roughly 

1 12 

the case that B eq oc pj . However, in all calculations below 
we compute B eq fully from Equation (0. In Section|4]we de- 
scribe observations that show the photospheric magnetic field 
strength B* is roughly linearly proportional to B eq for many 
stars. The measurements determine the constant of propor- 
tionality, and we use 



B* = 1135, 



eq 



(6) 

in the remainder of this paper. Figure |TJb) shows B* as a 
function of T e s and log g. 

We consider MHD waves that are driven by turbulent con- 
vective motions in the stellar interior. Th e original mod- 
els of wave gene r ation from tu rbulence (e.g., Lighthilll ll952t 
iProudmanl Il952t ISteinl 1 1967b dealt mainly with acoustic 
waves in an unmagnetized medium. More recently, how- 
ever, it has been shown that when a stellar atmosphere 
is filled with magnetic flux tubes, the dominant carrier of 
wave energy should be transver se kink-mode oscillations 
(Musi elak & Ulmschneide r 2002a). When the magnetic flux 
tubes extend above the stellar surface and expand to fill the 
vol ume, the kink-mode waves become shear Alfven waves 
(see lCranmer & van Ballegooiienll2005l). 

We u tilize the model results of iMusielak & Ulmschneiderl 
(2002a) to estimate the flux of energy in kink/Alfven waves 
in stellar p hotospheres. For simplicity , we use d only the sim- 
ulations of Mus ielak & Ulmschneiderl (l2002al) with their stan- 
dard parameter choices: a mixing length parameter of a = 2 
and a constant magnetic field strength that is 0.85 times the 
equipartition field strength. Our analytic fit to the results 
shown in their Figure 8 is 



"A* 



£eff 
7b 



exp 



Teff 

"7b" 



25 



where the dependence on g = log g is given by 
F 



10 9 erg cm 2 s 1 
7b 



5.724 exp 



11.48 



5.624 + 0.6002| . 



1000 K 

£ = 6.774 + 0.5057| . 



(7) 



(8) 

(9) 

(10) 
by 



These fits are simila r in form to those given 
iFawzy & Cuntzl (1201 lh for longitudinal MHD waves. 
Figure |2] shows a comparison between the above fitting 
formula and the plotted results of Mus ielak & Ulmschneiderl 
(2002a) for log g = 3, 4, and 5. Th e behavior of Fa* for lower 
values of log g was not given by Musi elak & Ulmschneiderl 
(2002a), bu t similar results were foun d for a wider range of 
gravities bv lUlmschneider et al.l (1 19961) for acoustic waves. 
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FIG. 2. — Comparison between the Musielak & Ulmschneider (2002a) nu- 
merical models (dot-dashed curves) and analytic fits (solid curves) for photo- 
spheric transverse wave energy fluxes Fa* as a function of effective tempera- 
ture and photospheric gravity. Numerical labels denote log g for each curve. 

We used the kink-mode energy flux to determine the trans- 
verse velocity amplitude vj_ of Alfven waves in the photo- 
sphere. The flux is defined as 



F A * = P*Vj_*Va 



(11) 



with Va* = B*/(47rp*) 1//2 being the photospheric Alfven 
speed. The above expression is not exact for waves under - 
going strong reflection (see, e.g jHeinemann & Olberj[l980l) . 
but it ends up giving a similar prediction for the height vari- 
ation of v± in the coron a that would come from a more ac - 
curate non-WKB model (ICranmer & van Bal legooijen 20051) . 
Figure 02c) shows how v^* varies as a function of T e ff and 
log g for solar metallicity stars. 

For the well-observed case of the Sun, we know that most 
of the photospheric magnetic field is concentrated into small 
(100-200 km diameter) flux tubes concent r ated in the in- 
tergranular downflow lanes (Solanki 1993; iBerger & Title! 
2001). The field strength in these tubes is close to equipar- 
tition, with B* « 1400 G. However, these flux tubes have a 
filling factor /* in the photosphere of about 0.1% to 1%, so 
the spatially a veraged magnetic flux den sity B*/* is only of 
order 1-10 G (ISchrijver & Harveyll989l) . 

2.2. Radial Evolution of Waves and Turbulence 

Figure [3] illustrates the stellar magnetic field geometry that 
we assume to exist above the surface of a cool star. Flux tubes 
that are open to the stellar wind 4 have a cross-sectional area 
A(r) that expands monotonically with increasing radial dis- 
tance r from the star. The condition V • B = demands that the 
product of A and the magnetic field strength B remains con- 
stant. Thus, B(r) inside a flux tube decreases monotonically, 
from its photospheric value of B*, with increasing distance. 
We normalize A such that at a given distance the total stellar 
surface area covered by open flux tubes is defined to be 



A = 4nr 2 f . 



(12) 



4 The presumed non-existence of magnetic monopoles implies that "open" 
field lines must eventually be closed far fro m the star, pr esumably via inter- 
actions with the larger-scale interstellar field i Davis 1955). 
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FIG. 3. — Summary illustration of flux tube expansion on a representative 
cool star. The filling factor (for open magnetic flux tubes) grows from /* in 
the photosphere to /tr at the transition region, and to an asymptotic value 
foo — ► 1 at large distances. The dimensions in this sketch are not drawn to 
scale. 

The dimensionless filling factor / tends to increase with 
height to an as ymptotic value of 1 as r — > oo (see also 
ICuntz et all 1 19991) . but its increase is not necessarily mono- 
tonic. We do not explicitly consider the properties of closed 
magnetic "loops" on the stellar surface, but the radial varia- 
tion of f(r) takes into account their presence. 

For each star, we intend to specify /» on the basis of ei- 
ther direct measurements or empirical scaling relations. The 
model described in Section 13.11 also requires specifying the 
value of / at the sharp transition region (TR) between the 
cool chromosphere and hot corona. We generally know that 
/* < /tr < 1, so in the absence of better information we will 
apply the assumption that that / TR = jf, where 6 is a dimen- 
si onless constant betwe en and 1. For the solar wind models 
of lCranmer et al.l (120071) the exponent 9 ranges between about 
0.3 and 0.5. 

Alfven waves propagate up from the stellar photosphere, 
partially reflect back down toward the Sun, de velop into 
strong MHD turbulence, and d issipate gradually (IVelli et al.l 
199 It IMatthaeus eTail 1 19991: ICranmer & v an Balleg ooiienl 



2005). Temporarily ignoring the reflection and turbulent cas- 
cade, the overall energy balance of an Alfven wave train is 
governed by the conservation of wave action. We define the 
flux of wave action S as 



/?Vj_Va(1 +M a ) 2 A = constant 



(13) 

where M\ = u /V\ is the Alfven Mach numb er and u is the 
radial outflow spee d of the wind (see, e.g., I Jacques! 119771; 
iTu & Marschl [T995). Close to the stellar surface, where 
Ma <C 1, this condition is equivalent to energy flux conser- 
vation (FaA = constant). In any case, the constant value of S 
in Equation ([TBI is known for each star because the conditions 
at the photosphere are known (and it is also valid to assume 
Ma — > there as well). The behavior of the wave amplitude 
as a function of density varies from vj_ oc p -1 / 4 close to the 
star to vj_ oc p +I / 4 at larger distances. 

The waves gradually lose energy due to turbulent dis- 
sipation, but for locations reasonably close to the stellar 
surface — e.g., the region shown in Figure [3] — it is not a bad 
approximation to use the undamped form of wave action 
conservation to compute the radi al dependence of vj_ (see 
ICranmer & van Ballegooiien 2005). Wave damping gives rise 



to plasma heating, and we adopt a phenomenological heating 
rate that is consistent with the total energy flux that cascades 
from large to small eddies. This rate is constrained by the 
properties of the Alfvenic fluctuations at the largest scales, 
and it does not specify the exact kinetic means of dissipation 
once the energy reaches the smallest scales. Dimensionally, 
it is similar to the rate of cas cading energy flux derived by 
Ivon Karman & Howarthl d 1 938b for isotropic hydrodynamic 
turbulence. The volumetric heating rate is given by 



Q = 



apv\ 



(14) 



(lHoflwegll 986: Hos sain et alTl 995: Zhou & Matt haeusH990t 
IMatthaeus et al.lll999h iDmitruk et~aT1l2002l) . The dimension- 
less efficiency factor a depends on the local degree of wave 
reflection and is discussed further below. The perpendicu- 
lar length scale is an effective correlation length for the 
largest eddies in the turbulent cascade. 

MHD turbulence occurs only when there exist counter- 
propagating Alfven wave packets along a flux tube. The star 
naturally creates upward waves, and we assume that linear re- 
flectio n gives rise to downward waves dFerraro & Plu mpton 
1958). We specify the ratio of downward to upward wave 
amplitudes by the effective reflection coefficient 1Z, and the 
efficiency factor a is given by 

^d+^)V2 

a = a ° (1+7^)3/2 < 15 > 

(see, e.g.. ICranmer et all 12007b . At the photospheric lower 
boundary, we assume total reflection with 1Z = 1 and thus 
a = ao- Higher in the stellar a tmosphere, we u se the low- 
frequency limiting expression of Cranmerl (1201 Ol) . 

Va-Moo 



K 



(16) 



Va + Moo 

where the wind's terminal speed is and we also assume 
that Ma <C 1 in the atmosphere. This expression also as- 
sumes that the wind speed at the point where Ma = 1 (pre- 
sumably far from the stellar surface) is roughly equal to u x . 
W e also set ap = 0.5 based on the turbulent transport models 
of iBreech et all d2009l) . 

3. MODELS FOR MASS LOSS 

In this section we present two complementary descriptions 
of cool-star mass loss that make use of the Alfven wave prop- 
erties discussed above. Supersonic winds can be driven by 
either gas pressure in a hot corona (Section lXTT i or wave pres- 
sure in a cool, extended chromosphere (Section [3~!2"l ). We first 
investigate each idea by assuming the other one is negligi- 
ble, and then we explore how to incorporate both processes 
together (Section l3~3l . 

3.1. Hot Coronal Mass Loss 

If the turbulent heating given by Equation ( TBI is sufficient 
to produce a hot (T > 10 6 K) corona, then the plasma's high 
gas pressure gradient may provide enough outward accelera- 
tion to produce a transition from a subsonic (bound) state near 
the star to a supersonic (outflowing) state at larger distances 
dParkerll 19581) . In this section we estimate the mass loss rate 
M of such a gas-pressure-driven stellar wind. 

We begin by computing 2* in the photosphere using p* and 
Vj_* in Equation ( TPfl i. For the Sun, we have the observational 
constraint that must be about the size of the granular mo- 
tions that jostle the flux tubes (i.e., roughly 100-1000 km). 
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For other stars we can assume that the horizontal scale of 
granulation r emains proportional to the photospheric pressure 
scale height (Robinson et al. 2004). Thus, we use 



A 



_L* 



Hq 



(17) 



where Hp, = 139 km and the models of 
ICranmer & van Ballegooiienl d2005l) were used to set 
the solar normalization of the correlation length to A^q = 300 
km. 

In the photosphere, we assume the turbulent heating is 
swamped by radiative gains and losses that are determined 
by the conditions of local thermodynamic equilibrium (LTE), 
and the temperature is set by those processes alone. At larger 
heights in the flux tube, the turbulent heating Q begins to have 
an effect. We define the chromosphere as the region in which 
Q is balanced by radiative losses. As one increases in height, 
however, the density drops to the point where radiative losses 
alone can no longer balance the imposed heating rate; this oc- 
curs at the sharp TR between chromosphere and corona. (See 
Section [32] for cases where this transition does not occur at 
all.) 

In the region between the photosphere and the TR, we as- 
sume that the wind flow speed is sufficiently sub-Alfvenic 
such that vj_ oc /T 1 / 4 . We also assume that A^ scales with the 
transverse size of the magnetic flux tube, so that Aj_ oc A 1 / 2 oc 
B~ l l 2 (Hollwe^mH). Thus, Equation (Hi can be rewritten 
as 

1/4 /Btr\ 1/2 

, hr (18) 



6tr 



Q!TR / Ptr 



where a* = 0.5 and all other photospheric quantities are as- 
sumed to be known. We also know that 



#tr _ f*_ ^ r\- 

rt ~ r ^ J* 



(19) 



B* /tr 

where the last approximation holds if there is a universal re 



lationship between /* and /tr as speculated in Section 12.21 
above. 

Just below the TR, the heating is just barely balanced by 
radiative cooling. In the optically thin limit, radiative cooling 
behaves as Q coo i = -n 2 A(T), where n is the number density in 
the fully ionized TR region. Let us then assume that Qtr = 
max|g C ooi|, where 



max|2 coo i| = 



PjR^ii 



(20) 



The quantity A max is the absolute maximum of the radiative 
loss curve A(T), and it occurs roughly at Ttr = 2 x 10 5 K. 
The value of A max depends on metallicity. To work out its 
dependence on Z, we computed a number of radiative loss 
curves for different metal abun dances using vers ion 4.2 of 
the CHIANTI atomic data base dYoung et al.ll2003h with col- 
lisional ionization balance (Mazzotta et al. 1998). We started 
with a traditional dGrevesse & Sauvall fT998) solar abundance 
mixture (Z/Z Q =1) and then recomputed A(T) by varying the 
metal abundance ratio Z/Z Q between and 10. We found 
that the maxima of the curves were fit well by the following 
parameterized function, 

7.4+42 • (21) 



10 23 erg cm 3 s 1 



Other examples of the metallicity dependence of A(T) 
have been given by, e.g., iBoehringer & Henslerl (1 19891) and 



iGnat & Sternbergl ( 12007!) . We have ignored any possible dif- 
ferences between a star's photospheric metal abundances and 
those in the low corona, although su ch differences have been 
measured in some cases (Testa 2010). 
With the above assumptions, we solve for the TR density, 



Ptr 



~ '/ 4 A 



4/7 



/ 



2(l-0)/7 



and we also derive the heating rate at the TR to be 



Q 



TR 



/ «tr6* 



8/7 



plA n 



1/7 



/ 



4(1-60/7 



(22) 



(23) 



A potential roadblock to solving Equations d22T - l23l is that we 
do not initially know the value of 5tr- This quantity depends 
on the reflection coefficient 1Z, which depends on the Alfven 
speed Va at the TR (see Equation (TToTl). which in turn depends 
on the unknown value of /?tr- In practice, we solve these 
equations iteratively. We start with an initial estimate of 1Z = 
0.5, we compute «tr, Ptr, and Va at the TR, and then we 
recompute 1Z for the next iteration. In all cases the process 
converges to a self-consistent set of values (with a relative 
accuracy of ~ 10~ 7 ) in no more than 20 iterations. 

The mass loss rate of the stellar wind is determined by the 
heating rate £>tr as well as other sources and sinks of energy 
at the TR. The general idea that the solar wind's mass flux 
is set by t he en e rgy balance at the TR w as first discussed by 
iHamm er ( 1982). lHansteen & Leerl d 19951) worked out the ba- 
sic scaling argument that is used below (see also iLeer et alj 
[T^IWthl^[T988t[Schwadron & McComasll2003l) . In the 
low corona and wind, the time-steady equation of internal en- 
ergy conservation is 



A dr 



FH-F coni +pu 



GM, 



. 



(24) 



where Fa is the energy flux associated with the heating, F con d 
is the energy flux transported by heat conduction along the 
field, and u is the outflow speed. The term in braces is con- 
stant as a function of radius, so it is straightforward to equate 
its value at the TR to its asymptotic value at r —> oo. The 
kinetic energy term proportional to u 2 is assumed to be negli- 
gibly small at the TR, but we assume it dominates the energy 
balance at large distances. Thus, 



-(puA) TR - 



GM* 



{puA)^ , (25) 



where F^jr is the heat flux Fh at the TR, and we realize 
that the product puA is also constant via mass flux conser- 
vation. We also make the key assumption that Uoo = V esc = 
(2GM*//?*) 1 / 2 , and thus we can write 



M 



puA = — — 2 (F Hi tr- 

esc 



(26) 



To evaluate Equation d26l l we need to estimate the value of 
fk.TR- Formally, Q = |V • Fh|, so to determine the magnitude 
Fh.tr one would have to integrate Q(r) along the flux tube. 
Taking account of the expanding flux tube area A oc Br 1 , and 
also assuming that r T R s» 



1 f°° 

Fh,tr = -— / drQ(r)A(r) 

A TR jR t 



(27) 
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Specifically, if Q oc r & and A oc r 1 , then 



|/3-7-l| 



QtrR*1i 



(28) 



cm 2 s 



Rather than specifying j3 and 7, we estimate the dimension- 
less scaling factor h by extracting both Q TR and F HT r from th e 
self-consistent solar wind models of ICranmer et al.l d2007l) . 
For a range of fast and slow solar wind solutions, we found 
that Qtr is typically between 1 .5 x 10~ 5 and 4 x 10~ 5 erg cm" 3 
and Fh.tr is typically between 8 x 10 5 and 3 x 10 6 erg 
This results in h usually being between 0.5 and 1.5. 
It is important to also verify that Fh,tr is less than the en- 
ergy flux carried "passively" by the Alfven waves as they 
propagate up from the photosphere. The latter quantity, which 
we call F a ,tr, represents the upper limit of available energy 
in the waves (at the TR) that can be extracted by the turbu- 
lent heating. Assuming that wave flux is conserved (i.e., that 
Ma < 1 at the TR), then F a , T r = /*F a *//tr. For the cool- 
star models discussed in Section [5] we found that the ratio 
Fh,tr/Fatr is usually around 0.1 to 0.5. Only in two cases 
did it exceed 1 (albeit with values no larger than 1 .5), and in 
those cases we capped Fh.tr to be equal to Fa.tr to maintain 
energy conservation. 

To evaluate Equation d26l i, we also need to estimate the 
magnitude of the downward conductive flux F con d from the ho t 
corona. For the solar TR and low corona, Withbroe (i 19881) 
found there to be an appr oximate balance between conduc- 
tion and radiation losses. Withbroe (1988) determined that 
F C ond ~ CmdPjR, where Ftr is the gas pressure at the TR, and 
the constant of proportionality is 



(29) 



where K e is the electron thermal conductivity, To s» 10 4 K is 
a representative chromospheric temperature, 7tr = 2 x 10 5 K, 
and c la d has units of speed. We evaluated the above integral to 
be able to scale out the metallicity-dependent factor given in 
Equation ( f2Tb above, and found that 



Crad 



14 J km,-' 



A max (Zn) 



(30) 



We used this expression to estimate F con( j = cwPtr. For the 
specific case of the Sun, conduction is relatively unimportant 
in open flux tubes, since F con( j < 0.05 Fh.tr- For the other stars 
modeled in this paper, the ratio F con d/FH.TR spanned several 
orders of magnitude from 10~ 4 to 10 _1 . In the eventuality 
that the estimated value of F con( j may exceed F h .tr, one would 
need an improved description of the coronal temperature T(r) 
to compute a more accurate value of the conduction flux. In 
our numerical code, however, we do not allow F con d to exceed 
a value of £Fh,tr, where £ is an arbitrary constant that we 
fixed to a value of 0.9. This condition was not met for any of 
the stars in the database of Section 

Once these energy fluxes are computed, we then compute 
M using Equations (1231 . d26l l. and d28l i. as well as the other 
definitions given above. Interestingly, this can be done with- 
out needing to know the temperature profile T(r). From a cer- 
tain perspective, the corona's thermal response to the heating 
rate Q may be considered to be just an intermediate step to- 
ward the "final" outcome of a kinetic-energy-dominated out- 



flow far from the star. 5 However, it should be possible to 
estimate the maximum coronal te mperature r m? x by invert- 
ing scaling laws given by, e.g ., lHammer et al.l d 19961) and 
ISchwadron & McComasl d2003l) . 

The mass loss rate given by Equation d26i > depends on our 
assumption that Uoo = V esc . Equation (l25T l shows that larger as- 
sumed values of would give rise to lower mass loss rates, 
and smaller values of Uoo would give larger mass loss rates. 
For the solar wind there is roughly a factor of three variation 
in Moo, from about 0.4V esc to 1 .3V esc . For other stars, it is rare 
to see observations where exceeds V esc , and Judge! d!992l) 
found generally that < V esc for luminous evolved stars. 
Even in the extreme case of = 0, however, Equation d25l 
would give only two times the mass loss assumed by Equa- 
tion d26l i. When compared to the larger typical observational 
uncertainties in M, factors of two are not a major concern. 

The Sun's mass loss rate of 2 x 10" 14 to3 x 1O" 14 M yr" 1 is 
modeled reasonably well with the model described here. The 
photospheric energy flux of Alfven waves is F A * w 1.5 x 10 8 
erg cm" 2 s" 1 , and the photospheric wave amplitude is vj_* ~ 
0.28 km s _1 dCranmer et alJl2007h . Although magnetogram 
observations sometimes give filling factors /* as large as 1 % 
at solar maximum, values of 0.1% tend to better represent 
the coronal holes that are conn ected to the largest volume of 
open flux tubes (see Figure 3 of Cranm er" & van Ballegooiienl 
120051) . Assuming /* = 0.001 and values for the other constants 
of ao = 0.5 and 6=1/3, Equations (I221423T ) give p T R ~ 5 x 
10~ 16 g cm" 3 and Qtr ~ 4 x 10~ 5 erg cm" 3 s" 1 , w hich are in 
agreement with the models of Cranmer et alJ d2007l) and oth- 
ers. Thus, with h = 0.5, Equation d26l i gives M sa 3.5 x 10" 14 
Af yr _1 . 

Although the above calculation of M is relatively straight- 
forward, it has not been boi l ed down to a simple scalin g 
la w such as that oflReimersI dl975l 1 19771) . iMullanl d!978l) . 
or lSchroder & Cuntzl d2005 ). However, if we make the fur- 
ther assumptions that a and h are fixed constants, and that 
Fh.tr 3> F con d, we can isolate several interesting scalings: 

1 . The ultimate driving of the wind co mes from the basal 
flux of Alfven wave energy F\*. Schro der & Cuntzl 
(2005) assumed that M scales linearly with F A *, but in 
our case we can combine the above equations with the 

12/7 

definition of 2* to find M oc F A J , which is noticeably 
steeper than a pure linear dependence. This positive 
feedback is qualitatively similar to what occurs in ra- 
diatively driven winds of more massive stars, for which 
the mass loss rate is proportional to the radiative flux (or 
luminosity) to a power larger tha n one (i.e., M oc Lj, 7 ; 
ICastorltanfT975l: IOwockill2004h . 

2. Extracting the dependence on magnetic filling factor, 

we found that M oc /* 4+3e) ^ 7 . Using the range of 9 from 
solar models (0.3-0.5), this gives a relatively narrow 
range of exponents, M oc 1 to f°- 8 . ISaarl d!996al) esti- 
mated that /* oc F r ~ J 8 for rotation periods F rot > 3 days 
(see also Section |4). Thus, for stars in the unsaturated 
part of the age-activity-rotation relationship, it may be 

5 Our approach, which ignores the details of this intermediate step, is an 
approximation that also sidesteps some other important issues. For example, 
a time-steady wind solution should pass through one or more critical points, 
and it and should also satisfy physical boundary conditions at r = R* and 
r — > 00. In Section|6]we summarize the necessary steps to producing more 
self-consistent versions of this model. 



THEORETICAL MASS LOSS RATES OF COOL STARS 



7 



£ 

o 

SB 
-. 



Of 



O 



-5 



•10 



•15 



hot corona 



giant 




cool chromosphere 



-25 



-5 



-20 -15 -10 

log P (g cm" 3 ) 

FIG. 4. — Density dependence of the heating rate Q defined in Equation 114) 
compared to th e de nsity dependence of the maximum radiative cooling rate 
from Equation {2Q\ (dashed curve). For the Sun, a numerical model (black 
solid curve) compares favorably with a simple analytic "bridging" between 
the near-star (Q oc p 1 / 2 ) and distant (Q oc p 9 / 4 ) scalings discussed in the text 
(black dotted curve). For an example evolved giant star (blue dot-dashed 
curve), the transition to the steeper density dependence occurs to the right of 
the cooling boundary. 



possible to estimate M oc f r 



■1.3 



These simple scaling relations are given only for illustrative 
purposes (see also Equation d45b below). The predictions of 
our "hot" coronal mass loss model should be considered to be 
the solutions of the full set of Equations ( fT7H30l l. 

3.2. Cold Wave-Driven Mass Loss 

In a high-density stellar atmosphere, it is possible that the 
turbulent heating described by Equation ( fT~4T > could be bal- 
anced by radiative cooling even very far from the star. In 
that case, hot coronal temper atures may never occur (see, e.g., 
ISuzukill2007t ICranmer 2008). The density dependence of the 
heating rate Q determines whether radiative cooling remains 
important at large distances, and Figure [4] shows two exam- 
ples of how Q may vary as a function of p. Near the stellar 
surface, where vj_ oc p" 1 / 4 and B oc p 1 / 2 , we can combine var- 
ious assumptions to estimate that Q oc p 1 / 2 . However, further 
from the star, where vj_ oc p +1//4 and B oc p, the density depen- 
dence becomes steeper, with Q oc p 9 ' 4 . 

Figure [4] compares the modeled heating rates with the 
"maximum cooling boundary" implied by Equation (1201 , The 
solar model crosses the boundary, and thus undergoes a tran- 
sition to a hot corona. The solid curve was taken from a nu- 
merical model of fas t solar wind from a polar coronal hole 
( ICranmer et alj 120071) . On the other hand, a model for a 
representative late-type giant sits to the right of the cooling 
boundary, which implies that radiative losses can maintain the 
circumstellar temperature at chromospheric values of ~ 10 4 
K even at large distances. The two models differ because 
the density at which Ma ~ 1 (where the Q curves undergo 
a change in slope) for the giant is several orders of magnitude 
larger than the corresponding density for the solar case. This 
density is an output of a given mass loss model and cannot be 
specified a priori. 

In this section we develop a model for cool-star mass loss 
under the assumption of strong radiative cooling. In this 



case the lParkerl(fl958l) gas pressure driving mechanism cannot 
drive a significant outflow. However, when the flux of Alfven 
waves is large, they can impart a strong bulk acceleration to 
the plasma due to wave pressure, which is a nondissipative net 
ponderomotive force exerted by virtue of wave propagation 
through an inhomogeneous medium ( Bret herton & G arrett 
Il968t lJacqueslll977l) . The subsequent calculation of M for 
a "cold wa ve-driven" s t ellar w ind larg ely follows the devel- 
opment of iHolzer et al.1 (1 19831) (see also ICranmerll2009l) . 

Three key assumptions are: (1) that the Alfven wave ampli- 
tudes in the wind are larger than the local sound speeds, (2) 
that there is negligible wave damping between the stellar sur- 
face and the wave-modified critical point of the flow, and (3) 
that the critical point occurs far enough from the stellar sur- 
face that / ~ 1 there (i.e., the flu x tube expans i on bec omes ra- 
dial). A fourth assumption from lHolzer et al.1 (119831) — which 
was initially not applied here but later found to be valid — 
is that the stellar wind is sub-Alfvenic at the critical point 
(i.e., that M A <C 1 at the critical point). ICranmeri d2009l) con- 
structed a set of numerical models for the cold polar outflows 
of T Tauri stars that did not make the fourth assumption, and 
found that for a wide range of parameters the assumption was 
justified. Thus, here we u se the value of th e critical radius 
given by Equation (35) of IHolzer et al.l (119831) . in which all 
four of the above assumptions were applied, and 



'"crit 



7/4 



l+(vWKsc) 2 



(31) 



Once the critical point radius is known, it becomes possible 
to use the known properties of the Alfven waves t o determine 
the win d velocity and density at the critical point. Holz er et~aT] 
( 1983) found analytic solutions for these quantities in the lim- 
iting case of Ma <C 1 at the critical point. Here we describe 
a slightly more self-consistent way of computing Ma and the 
mass loss rate, but we also continue to use Equation OTb that 
was derived in the limit of Ma <C 1 . There are three unknown 
quantities and three equations to constrain them. The three 
unknowns are the critical point values of the wind speed u, 
density p, and wave amplitude vj_. The first equation is the 
constraint that the right-hand side of the time-steady momen- 
tum e quation must sum to zero at the critical point of the flow 
(e.g.. iParkerll 19581) . For the conditions described above, this 
gives 



2ui 



GM* 



. 



(32) 



and it is solved straightforwardly for u crlt . The second and 
third equations are, respectively, the definition of the critical 
point velocity in the "cold" limit of zero gas pressure, 



1 + 3M A 
1+M A 



(33) 



and the conservation of wave action as given by Equation 
( fT3l . The fact that Va appears in these equations and depends 
on the (still unknown) density makes it difficult to find an ex- 
plicit analytic solution for p CI j t . We again use iteration from 
an initial guess to reach a self-consistent solution for u, p, and 
Vj_ at the critical point. The stellar wind's mass loss rate is 
thus determined from M = 47rr 2 rit M cl -j t /9 C i-it- 

Because the mass loss rate is set at the critical point, we 
do not need to specify the terminal speed Koq. For most im- 
plementations of the above model, the denominator in Equa- 
tion OTI ) is close to unity and thus we have M 2 rit V^ c /7, 
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or that M cl -i t is about 38% of the presumed value of u^. 
Of course, there have been stellar wind models with non- 
monotonic radial variations o f u(r), with M crrt > Mqo (e.g., 
iFalceta-Gon calves et al. 2006). It is also possible for "too 
much" mass to be driven past the critical point, such that 
parcels of gas may be decelerated to stagnation at some height 
above r cl -i t and thus would want to fall back down towards the 
star. In reality, this parcel would collide with other parcels that 
are still accelerating, and a stochastic collection of shocked 
clumps is likely to result. Interactions between these parcels 
may result in an extra degree of collisional heating that could 
act as an extended source of gas pressure to help maintain a 
mean net outward flow. Situations similar to this have been 
suggested to occur in the outflow s of pulsating cool stars 
dBowenl Il988t IStruck et al.1 120041) . T Tauri stars (j Cranmer 
120081) . and luminous blue variables ( Ivan M arie et al. 2009^ 

3.3. Combining Hot and Cold Models 

A proper treatment of a stellar wind powered by MHD 
turbulence — and accelerated by a combination of gas pressure 
and wave pressure effects — requires a self-con sistent numeri 



cal solution to the cons ervation equa t ions (e . g.JCranmer et al 
20071; ISuzukil l2007t ICohen et al.1 120091; lAirapetian et al 



20101) . However, in this paper, we explore simpler ways of 
estimating the c om bine d effects of both processes. 

Sections [3.11 and l3~2l gave us independent estimates for the 
mass loss rate assuming only gas pressure or wave pressure 
were active in the flux tube of interest. We refer to these two 
mass loss rates as Mhot and M co id, respectively. It seems clear 
that when one of these values is much larger than the other, 
then one process is dominant and the actual mass loss rate 
should be close to that larger value. For the manifestly "hot" 
example of the Sun, we found that Mhot/M co id ss 20, which 
correctly implies that gas pressure driving is dominant. For 
most examples of late-type giants with L„ > 10L Q , the ratio 
Mh ot /M co \i was found to decrease to values between about 0. 1 
and 3. This could mean that gas and wave pressure gradients 
are of the same order of magnitude for these stars. 

One of the most straightforward things that can be done 
is to assume the combined effect of gas and wave pressure 
produces a mass loss rate equal to the sum of the two in- 
dividual components, Mhot + M co id. This preserves the idea 
that one dominant mechanism should determine M when the 
other would predict a negligibly small effect. It also makes 
sense based on Equation d24i i. which shows how the energy 
fluxes sum together linearly in the internal energy equation. 
If there were multiple sources of input energy flux, Equation 
d26l i would show that the resulting mass loss rate should be 
proportional to their sum. 

However, there is one complication that hinders us from 
simply adding together Mhot and M C oid- The calculation of 
Mhot from Section 13.11 contains the assumption that the TR 
turbulent heating always obeys the near-star density scaling 
Q oc p l l 2 . It therefore predicts that all stars eventually undergo 
a transition to a hot corona. For some stars (like the late-type 
giant in Figure 0), however, we know that there should be no 
corona and it is erroneous to assume that Mhot has any real 
meaning. Thus, for each model we compute the wind speed 
at the TR from mass flux conservation, 

M ho 



tion, the TR Mach number Ma.tr = wtr/Va,tr should be much 
smaller than one. As expected, this condition was found to 
be violated for late-type giants having L* > IOOLq. Thus, in 
these cases we should replace Mhot with either a drastically re- 
duced value or zero — the latter in cases where the Q{p) curve 
always falls to the right of the maximum cooling boundary 
in Figure [4] After some experimentation, we found that re- 
ducing the initially computed value of Qtr by a factor of 
exp(-4M ATR ) does a reasonably good job of reproducing the 
result of using a more consistent Q(p) function. Thus, we pro- 
pose that the summing of the "hot" and "cold" mass loss rates 
be done with the following approximate expression, 

M ss M C oid+Mhotexp (-4M AXR ) (35) 

where Mhot and Ma.tr are computed using the assumptions of 
Section [3~T1 and M co id is computed using the model given in 
Sectionl3~2l 



4. MAGNETIC ACTIVITY AND ROTATION 

An important ingredient in the above models — which re- 
mains unspecified for most stars — is the photospheric filling 
factor /*. It is now well-known that both /* and the mag- 
netic flux density B* f * exhibit significant correlations with 
stellar rotation speed dSaar & Linsky|Jl 986; Ma rcv & Basril 
1989; Montesinos & Jordanll993l:ISaarl2001l) . For many stars 
the rotation rate also scales w ith age, chromospheric activity, 
and coronal X-ray emission (Skumani chl 119721; iNoves et alj 
[l984t IPizzoiato et al]l2003t iMamaiek & Hillenbrand! 120081) . 
A prevalent explanation for these correlations is that an 
MHD dynamo amplifies the magnetic flux in proportion 
to th e large-scale energy input from differential rotation 
(e.g., |Parkgi [192 91; Montesin os et al.1 l200 1 j iBushbvl 120031: 

niiikj - 



20091: HsiketaD 



MTR = zw 1 — ' (34) 

and we demand that for Mhot to have a consistent interpreta- 



Mos s & Sokolol I2009t IChristensen et all 
120111) 

In this section we construct an empirical scaling relation 
that will allow a reasonably accurate determination of /* as 
a function of P m and the other stellar parameters. Other 
estimates of this relationship have been made in the past 
(iMontesinos & Jordanl 19"93HStepienl 1994l:ISaaill996aLl2001h 
ICuntz et all Il998t iFawzy et alJ 120021) . However, since our 
goal is to apply this relation to stellar wind acceleration (in 
open flux tubes that cover a subset of the inferred /* area) and 
to evolved giants (which are greatly undersampled in obser- 
vational studies of /*), we aim to reanalyze the existing data 
rather than rely on other published scalings. 

Table 1 lists the properties of 29 stars that have reliable 
measurements of their fundamental parameters, rotation rates, 
and either independent or combined values of B* and /*. The 
sources for these values are given as numbered references 
in the final column. In many cases the available sources 
gave only a subset of the basic stellar parameters. When 
necessary, we used Equation (Q]i and information from the 
NASA/IPAC/NExScI Star an d Exoplanet Database (NStED) 6 
to fill in missing values (see Berriman et al.ll20loh . Table 1 
also gives approximate "quality factors" q that describe the 
relative accuracy of the measurements, and in the Appendix 
we describe these factors in more detail. 

It has been known for some time that, for dwarf stars, 
B* never appears to be very far from t he equipartition field 
strength B eq (e.g.. ISaar & Linsky| [l986). Figure [5] plots the 
ratio B*/B eq for the measurements in Table 1 that have sepa- 
rate determinations of B* and /* . The sizes of the symbols are 

6 |http : / /nsted. ipac . caltech . edu7] 
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Table 1. Magnetic Activity Data for G/K/M Stars 



Name 


Trff (K) 


log £ 


M* /Mq 




L*/L Q 


Prot (d) 


[Fe/H] 


B* x /* (G) 


Qual. 


Ref. 


Sun 


5770 


4.44 


] 


j 


\ 


25.3 


o 


1 4D0 y (0 00 1 -0 o n 






59 Vir (GO Vl 


6234 


4.60 


1.17 


0.897 


1. 10 


3.3 


+0.280 


1000 y 19 

1 \.J\J\.J A \J . 1 y 


2 


12 3 


















500 


4 


1, 2, 3 


■v 1 Ori (GO V> 

Y Wll \ VJU V / 


5955 


4.30 


0.67 


0.962 


1 Q5 


5.2 


-0.039 


1 000 y 60 


1 


4 5 6 7 


1J >Jji>- 1 vJ 1 V f 


5905 


4.30 


0.8 1 


1 .05 


1 .21 


1 3.5 


+0.024 


1 800 y 1 
i cj\j\j a \j . i w 


1 


8 5 7 


5R Fri fni 

JO LI 1 i 1 V J 


5826 


4.54 


1.02 


0.898 


0.838 


10.8 


-0.013 


330 


3 


9 10 7 


q rw ("n? v i 

y y^ c i l uz v ^ 


5790 


4 40 


1 


1 04 


1 1 1 


7 7 


+0 159 


1400 y 39 

1 H\J\1 A U.Ji 


1 


4 5 6 7 


un 9soqq /"n? v». 


5761 


4 37 


94 


1 05 


I 10 


8 7 


+0 137 


1 700 y 30 


1 


115 6 7 


k- rw (ns vi 

ri* v^-Cl ^VJJ V ) 


5771 


4.56 


1.02 


0.877 


0.770 


9.4 


+0.056 


321 


2 


9 5 7 


















392 


2 


9, 5, 7 


















406 


2 


9 5 7 


















480 


2 


9 5 7 


















1 500 y "15 

1 -J\J\J A U. JJ 


\ 


9 5 7 


DUU rt ( VJ O V ) 


5551 


4 57 


86 


801 


550 


6 2 


—0 122 


1 600 y 99 

1 UVAJ A XJ.Z.Z 


2 


12 5 7 


















1 £00 y "IS 

1 OVA./ A U.JJ 


1 


12 5 7 


















9000 y 70 


2 


12 5 7 


















1900 y 18 

1 j\j\.J A U.lu 


3 


12 13 7 


un 1 S93Q 1 rtTfl. S V*l 


5495 


4 30 


86 


1 09 


971 


1 1 1 


—0 049 


1 700 y 1 S 

1 / A U. lo 


1 


4 5 6 7 


/ \J \jVjvi r\. LAV V ^ 


5300 


4.52 


0.89 


0.86 


0.53 


19.7 


+0.040 


1900 y 18 

1 Lajki A U.lu 


1 


14 5 


DF Ron (PC 1 \H 


5231 


4.45 


0.78 


0.871 


0.512 


9,0 


+0. 108 


1 700 y 06 

1 l\j\.J A \J.\tYJ 


2 


3 6 7 


HP) 1 7975 ("K? 
sni_j i / j \ iv^ v ^ 


5225 


4.40 


0.79 


0.93 


0.58 


6,6 


+0.067 


1 500 y 1 5 

1 -J\J\.J A U.JJ 


3 


4 5 6 7 


LQ Hya (K2 V) 


5070 


4.68 


0.80 


0.673 


0.270 


1 .6 


+0.330 


3500 x 0.70 


2 


15, 3, 7 


un i i S404 fK? 


4814 


4.53 


0.77 


0.791 


0.303 


1 8.8 


-0. 193 


9 1 00 y 90 


1 


16 5 6 7 


HD 4628 (K2.5 V) 


5004 


4.64 


0.77 


0.69 


0.273 


38.5 


-0.27 


1600 y \1 


] 


17, 6, 7 


OTT Gem (K^p Vl 

\_f LJ vJClll 1^I\JC V f 


4959 


4.30 


0.610 


0.915 


0.457 


7.4 


-0. 170 


9400 y 50 

Z.t\J\./ A U.Ju 


\ 


8 5 7 


fc JELll IJV+.J V ) 


5094 


4 60 


83 


754 


345 


1 1 7 


—0 097 


165 


3 


18 5 10 7 


















1 000 y ^0 


1 


18 5 10 7 


















1900 y \1 


2 


18 5 10 7 


















1 440 y OSS 

I'rtU A U.WOO 


4 


18 13 10 7 


36 Dnh A (K5 W\ 


5 135 


4.54 


0.602 


0.69 


0.299 


20.3 


-0.206 


1 500 y 13 

1 JUU A U. 1J 


1 


19 5 10 7 


3fS Onh R (V> V\ 

JU WU11 JE> (IV J V ^ 


5 103 


4 58 


486 


59 


213 


22 9 


—0 195 


60 


2 


19 5 10 7 


VS33 Tan (K"Sp Vl 

V OJJ Id. LI ^IvJC V ) 


4450 


4 57 


80 


77 


209 


1 85 


+0 340 


9600 y SO 


3 


20 21 5 7 
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proportional to the observational quality factors, and all statis- 
tical fits and moments discussed below were weighted linearly 
with q. Figure 13a) shows that th ere is no strong correlation 
of B*/B eq with r e ff. Saar (Il996ah found a slight increase in 
B*/B eq for the most rapid rotators (P mt < 3 days), and Figure 
13b) shows that when more data are included this trend sur- 
vives but is not strong. The power-law fit is consistent with a 
relationship B* oc 3 , but we do not consider it significant 
enough to apply it below or to extrapolate it to longer rotation 
periods. 

We found that the ^-weighted mean value of Z?„/B eq for the 
entire sample (1.16, with standard deviation ±0.38) is only 
marginally higher than the mean value for the subset of slower 
rotating, non-saturated stars with P mt > 3 days (1.13, with 
standard deviation ±0.25). Equation (O gives the latter mean 
value, which we use in Section [5] for modeling the winds of 



the (generally slowly rotating) stars with observed mass loss 
rates. We also use Equation © to estimate /* for the cases 
where only the product B*/* has been measured. 

A primary indicator of stellar magnetic activity appears to 
be the photospheric filling factor . There have been a num- 
ber of different proposed ways to express the g eneral anti- 
correla tion between activity and rotation period. Noy es et aD 
(1984) found that indices of chromospheric activity corre- 
late better with the so-called Rossby number Ro = P mt /T c , 
where t c is a measure of the convective turnover time, than 
with P rot alone. For other data sets, however, the useful- 
ness of the Rossby number h as be e n call ed into question 
dBasril 119861: IStepienl 11994 . ISaari (1199 ll) postulated that 
fi*/* (an d presumably also f t itself) is proportional to Ro" 1 
(see also iMontesinos & Jo rdan 19931 ICuntz et alJll998l : ls"aarl 
l200ltlFawzvetal.ll2002l) . 
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data (red) and for only stars having P rot > 3 days (blue). Dotted and dashed 
lines show regions within ±1<t of the means. Quality factors are denoted by 
crosses (q = 1), triangles (q = 2), squares (q = 3) and circles (q = 4). 

To compute the Rossby number for a given star, we need 
to know the convective turnover time r c . Figure [6] compares 
several past calculations of r c with one another. For most stars 
we will utilize a parameterized fit to the set of ZAMS stellar 
models given bv lGunn et ail (fl998l) . 



t c = 314.24exp 



'eff 



1952.5K 



'eff 



6250 K 



+ 0.002, 
(36) 

where t c is expressed in units of days and the fit is valid for the 
approximate range 3300 < r e ff < 7000 K. Such a fit ignores 
how t c may depend on other stellar parameters b esides effec- 
tive temperature, but more recent sets of mode ls (Landi net al] 
I2010t IBarnes & Kiml |2010t iKitchatinov & Olemskovll201 lb 
also found reasonably monotonic behavior as a function of 
r e f t for a broad range of stellar ages and masses. 

There are indications that the simple relationship between 
T c and T e s seen for main-sequence stars is not universal. For 
example, 

1. Low-mass M dwarfs (with < 0.35 M ) are likely 
to be fully convective, and thus their dynamos are 
likely to be driven by fundament ally different processes 
than exist in more massive st ars (Mullan & Mac Donaldl 
120011: iReiners & Basrill2007t Hrwin et alj|201 lb . There 
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FIG. 6. — Estim ates of the convecti ve turnover time t c as a fun ction of T e ff. 
ZAMS model s oflGunn et all JT998h were fit by Equation {36) (black solid 
curve). The IGunn et al] 0.998) evolutionary track for a 2.2Mq star (blue 
dashed curve) outlin es an upper limit for r c at intermediate temperat ures. Lo- 
cal val ues of t c from Landin et al. 1 2010) (green symbols) and Barn es & Kiml 
1 2010) (orange dot-dashed curve), the parameterization given by N oyes et aH 
1 1984) (black dotted curve), and the M dwarf estimate of t c « 70 d used by 
Reiners et al. (2009) (red box labeled by "dM") are also shown. 

is also some disagreement about the relevant t c val- 
ues for these s t ars. F igure [6] shows that the models of 
IBarnes & Kiml (12010b exhibit a sl ight discontinuity a t 
the fully convective boundary. The Rei ners et al.l (2009) 
semi-empirical estimate of t c m 70 days for M dwarfs is 
about a factor of 2-3 lower than t hat of | Barnes & Ki 
d2010b . However, because the IReiners et al.1 d200' 
value is in reasonable agreement with an extrapolation 
of Equation (f36t to lower effective temperatures, we 
will just use this expression and not make any special 
adjustments to the Rossby numbers of fully convective 
M dwarfs. 

2. Luminous evolved giants exhibit qualitatively different 
interior properties than do main sequence stars of sim- 
ilar r e ff. Despite not having firm measurements of the 
magnetic activities of evolved giants, we will want to 
estimate /» for such s t ars in order t o compute their mass 
loss rates. iGondoinl d2005l l2007b found that the cor- 
relation between X-ray activity and rotation in G and 
K giants is consistent with that of main-sequence stars 
if the larger valu es of r c from the evolved models of 
iGumTetaD (fl998h were used instead of the ZAMS val- 
ues! (seejhe blue dashed curve in Figure |6]). Similarly, 
lHalll (1 1994b calculated luminosity-dependent scaling 
factors that can be used to multiply the ZAMS value 
of t c to obtain a consistent relati on between ro tation 
and photometric activity (see also lChoi et alJl995l) . We 
found that the above results can be generally repro- 
duced by multiplying the ZAMS value of r c by a fac- 
tor (go/g) ' 23 , which applies only for low-gravity sub- 
giants and giants (i.e., only for g < g@). In Section [5] 
we explore the extent to which this kind of approximate 
correction factor helps to explain the activity and mass 
loss of evolved stars. 

A slightly different way of estimating the magnetic flux of 
a rotating star is to take advantage of a proposed "magnetic 



THEORETICAL MASS LOSS RATES OF COOL STARS 



11 



Bode's law;" i.e., the conjecture that the star's ma gnetic mo- 
ment scales linearly with i ts angular momentum ( Ar ge et all 
119951; [Baliunas et al.ll 1996b - Using the stellar parameters de- 
fined above, this corresponds approximately to 



B*f*Rl cx 



MJLi 



(37) 



The above relationship does not take into account variations 
of the moment of inertia (for different stars) away from an 
idealized scaling of / ~ M*Rl, and it assumes the magnetic 
moment is dominated by a large-scale dipole component. It 
is possible to test this idea with the data given in Table 1 by 
evaluating the correlation between /* and M^iR^P^B^y 1 . 

Figure [7j shows how the empirical set of /» values corre- 
lates with rotation period, Rossby number, and the proposed 
magnetic Bode's law. Rather than use the Rossby number it- 
self, we instead plot the data in Figure 13b ) as a function of a 
normalized ratio Ro/Ro©, where according to Equation d361 l. 
the Sun's Rossby number Ro Q = 1.96. Such a normalization 
allows us to neglect any scaling discre pancies between "lo- 
cal" and "global" d efinitions of r c (e.g.. lPizzolato et afll2001l: 
lLandin et al.l l2010b . The Sun's large range of measured /* 
values (10~ 3 to 10~ 2 ) is indicated with a vertical bar, and it 
is likely that all other stars exhibit such a range on both rota- 
tional and dynamo-cycle time scales. 

Figures [71 a) and |7Jc) show that the correlations with P rot 
and the proposed magnetic Bode's law are not especially 
strong. However, if all of the lowest quality (q = 1) mea- 
surements were removed, the correlation with P rot would be 
improved significantly. Figure [71b) shows that the Rossby 
number seems to be a slightly better ordering parameter, and 
it compares the individual data points with several functional 
relationships. The blue and red solid curves are subjective fits 
to the minimum and maximum bounds on the envelope of data 
points, with 

0.5 

/min — F1 . ; lri ^xi^i -x 7 (38) 



fmax ~ 



[l + (x/0.16) 26 ] 13 ' 
1 



l+(x/0.31) 2 - 5 



(39) 



where x = Ro/Roq. W e also show empirical and th eoreti- 
cal fitting formulae from Montesinos & Jordan ( 1993). Other 
comparison s could also be made with r elationships given by 
ICuntz et al.1 (119981) . iFawzy et all (120021) . and others, but they 
all appear to fall near the green and red curves. 

Note that for slow rotation rates (i.e., large Rossby num- 
bers) the scaling laws shown in Figure |7jb) imply a sig- 
nificantly steeper decline of /* than has been suggested in 
the past. For e xample, Saari (119911) estimated /» cx Ro -1 , 
and lSaarf d!996al) estimated /* cx P^ t & . On the other hand, 
our empirical upper and lower bounds suggest / max cx Ro~ 2 5 
and /mji, cx Ro~ 3 4 respectively. This is similar to the ob- 
serv ed relationship between R ossby number and X-ray activ- 
ity. iMamaj ek & Hil lenbrand! (120081) found that the ratio of 
X-ray to bolometric luminosity Lx/Lboi drops by about a fac- 
tor of 700 as the Rossby num ber increases by a factor of ten 
from 0.25 to 2.5 (see also IWright et alJl201 ll) . This corre- 
sponds very roughly to a power-law decrease of Ro~ 2 85 . Its 
agreement with the behavior of /* shown above is also consis- 
tent with existing empirical corr elations between X-rays and 
magnetic activity dPevtsov et al.ll2003h . 

In addition to the rotational scaling of /* with Ro, 
there is also likely to be a "basal" lower limit on the 
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FIG. 7. — Comparison of possible correlations between measured /* rilling 
factors with (a) rotation rate, (b) Rossby number, and (c) a magnetic Bode's 
law parameter (see Equation 137) ). Solid curves in (b) denote the lower (blue) 
and upper (red) envelopes surrounding the data, and green curves show fitting 
fo rmulae from Equations 2.3 (do tted; empirical) and 7.3 (dashed; theoretical) 
of Montesinos & Jordan ( 1993). Quality factors are denoted by the same 



symbols used in Figure[' 
bar. 



and the Sun's range of /» is shown with a vertical 



outer atmo spheric activity of a sta r (e.g.. ISchriiverl Il987h 
Cuntz et alJll999tlBercik et al.ll2005t l Takeda & Takada-Hidail 
201 U iPerez Martinez et all 1201 lb . Whether this lower limit 



is the result of acoustic waves, a turbulent dynamo, or some 
other physical process, there is probably a minimum value 
of /» that is inde pendent of rotation rate. For example, 
iBercik et al.l ([2005) found that turbulent dynamos in main se- 
quence stars can generate a flux density fi*/* k4G without 
much variation from spectral types F0 to M0. Using Equa- 
tion (O for fi*, we can estimate a basal filling factor for these 
stars of /» 0.001-0.002, which is close to the Sun's min- 
imum value. However, since it is still uncertain whether or 
not the Sun has exhibited truly basal flux conditions in recent 
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years deliver & LindlloTlh . we will set a slightly lower value 
of /basal = 10~ 4 to be used in the mass loss models below. 

Before moving on to apply the empirical values of /* to our 
model of mass loss, we emphasize that the measurements do 
not directly provide the filling factor of open magnetic flux 
tubes. Ideally, Zeeman broadening measurements should be 
sensitive to the total flux in strong magnetic elements on the 
stellar surface, no matter whether the field lines are closed 
or open. In many cases, however, the closed-loop active re- 
gions have significantly stronger local field strengths than the 
open regions. Therefore the closed-field regions are likely to 
dominate the spectral line broadening that gives rise to the 
observational determinations of /» (see the Appendix). With- 
out spatially resolved magnetic field measurements, we do not 
yet have a definitive way to predict how a g iven star divides 
up its flux tubes between open and closed. Mes tel & Spruitl 
( 1987) claimed that as the rotation rate increases (from slow 
values similar to the Sun's), the relative fraction of closed 
field regions should first increase, then eventually it should 
decrease as centrifugal forces strip the field lines open. We 
can speculate that the spread in the measured /* data may tell 
us something about the closed and open fractions. Because 
closed-loop active regions tend to have stronger fields than 
open coronal holes, the lower and upper envelopes that sur- 
round the data in Figure [7Jb) could be good proxies for the 
filling factors of open and closed regions, respectively. More 
specifically, we hypothesize that f m [ n is seen when no active 
regions are present on the visible surface (i.e., / m ; n s» / ope n) 
and that / max is seen when active regions dominate the ob- 
served magn etic flux. This idea is tested, in a limited way, in 
Section EJ 

5. RESULTS 

Here we present the results of solving the mass loss equa- 
tions derived in Section [3] using the empirical estimates for 
the rotational dependence of the magnetic filling factor de- 
rived in Section @] For hot coronal mass loss, we assumed 
values for the dimensionless parameters ao = 0.5, h = 0.5, and 
9=1/3. As discussed above, these values were "calibrated" 
from our more detailed knowledge of the Sun's coronal heat- 
ing and wind acceleration. Our use of these values for other 
stars is an extrapolation that can be tested by comparison with 
observed mass loss rates. 

5.1. Database of Stellar Mass Loss Rates 

Figure |8]is a broad overview of observed stellar mass loss. 
It plots the locations of individual stars in a Hertzsprung- 
Russell type diagr am with their mass lo ss rates shown as sym- 
bol color (see also lde Jager et al.lll988l) . A box illustrates the 
approximate regime of parameter space covered by the mod- 
els developed in this paper; it extends from the main sequence 
up thr ough the regime of the so-called "hybrid chromosphere" 
stars (lHartman n et alj[l980l) . and possibly also into the pa- 
rameter space of cool luminous supergiants. In addition to the 
cool-star data discussed below, we also includ e in Figure HI 
measured mass loss rates of hot, massive stars dWaters et all 
1987 1 ; lLamers et al.1 119991; Mokiem et all ]2007t ISearle et all 
20081). FGK supergiants (ide Jager etaijll988h . AGB stars 
(Berge at & Chev allier 2005; Guandalini 120 1 Oh . red giants in 
globular clusters dMeszaros et al.l2009l). and M dwarfs in pre- 
cataclysmic variable binaries (De beil2006l) . Many of these 
stars are not included in the subsequent analysis because we 
have no firm rotation periods or magnetic activity indices for 
them. 



Table 2 lists the properties of 47 stars for which our knowl- 
edge appears to be complete enough to be able to compare 
theoretical and observed values of M. For the Su n, the 
range of volume-integrated mass loss rates comes from lWangl 
(1998). The sources for all listed values are given as num- 
bered references that continue the sequence started in Table 
1; the citations corresponding to numbers 1-32 are given in 
Table 1 . In cases where T e a, log or [Fe/Hl were estimated 
from the PASTEL database (Soub iran et alJl2010l) . we aver- 
aged together multiple measurements when more than one 
was given. In the few cases where the same star appears in 
both Table 1 and Table 2, for consistency's sake we will re- 
compute /* from the star's rotation p eriod when calculating 
theoretical mass loss rates (see Section l5T2l ). 

For binar y systems _ with astrospheric measurements of M 
(see, e.g., IWood et alll2002h . the numbers given are assumed 
to be the sum of both stars' mass loss rates. We list that same 
value for both components and denote it with "(A+B)." We did 
not utilize the published astrospheric measurements of Prox- 
ima Cen and 40 Eri A, which gave only upper limits on M, 
and A And and DK UMa, which had uncertain detec tions of 
astrospheric H I Lya absorption dWood et alj|2005alb1) . 

At the bottom of Table 2 we list three stars that have pa- 
rameters at the outer bounds of what we intend to model. 
They are test cases for the limits of applicability of the 
physical processes summarized in Section [3] EV Lac is 
an active M dwarf and fl are star that proba bly has a fully 
convective interior (e.g., lOsten et all 120101) . Such stars 
may exhibit qualitatively different mechanisms of mass loss 
and rotation-activity correla t ion than do stars higher up 
the main sequence dMullanl fl996t iReiners & Basil 120071: 
Irwin etalJl20lU iMartmez-Arnaiz et al.1 1201 lb IVidotto et ail 



201 11) . V Hya is an N-type carbon star with an extended 



and as ymmetric AGB envelope and evidence fo r rapid ro- 
tation dBarnbaum et alJ 11995b iKnapp et al.1 Il999t) . 89 Her 
is a post-AGB yellow supergiant with multiple detections 
of circumstellar neb ular material (Sargent & Osmerl Il969t 
Bujarrabal et al. 2007). It is worthwhile to investigate to what 
extent the mass loss mechanisms proposed in this paper could 
be applicable to these kinds of stars. 

Not all stars in Table 2 have precise measurements for their 
rotation period. For 61 Vir and 70 Oph B, we used published 
estimates of the rotation period that were obtained from the 
known c orrelation between ro tation and chromospheric Ca II 
activity (Baliun as et alJfT996l) . For essentially all stars more 
luminous tha n ~ 5L Q (with the exception of HR 6902; see 
GriffiniCijll) we estimated P mt via spectroscopic determina- 
tions of vsin/ from the rotational broadening of photospheric 
absorption lines. The inclination angle / is the main unknown 
quantity. Chandrasekhar & Munch d!950h found that for an 
isotropically distributed set of inclination vectors, the mean 
value of sin / is ir/4. Thus, we estimate a mean rotation pe- 
riod 

<^> = 717^- • (40) 
(4/7r) vsin; 

Note that the median of sin/ for an isotropic distribution is 
not equal to the mean; the former is given by In order 

to encompass both values, as well as the majority of "most 
likely" values of f rot , we can adopt generous uncertainty lim- 
its for which we will estimate /* and the other derived quan- 
tities for mass loss. For the isotropic distribution of direc- 
tion vectors, the quantity sin i falls between 0.5 and 1 approx- 
imately 87% of the time. This is a reasonably good definition 
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FIG. 8. — Hertzsprang-Russell diagram showing observed mass loss rates. The color scale at lower-left specifies log M, where M is measured in Mq yr . An 
estimate for the ZAMS is shown in gray, and the approximate domain of parameter space covered by the models of this paper is outlined by a black dotted box. 



for uncertainty bounds that would correspond to ±1 .5cr if the 
distribution were Gaussian. Thus, for stars with only vskn 
measurements, we use the following values as error bars on 
the derived rotation period: 



2 ^ P m t <- 4 



(41) 



5.2. Comparing Predictions with Observations 

We applied the combined model for mass loss that culmi- 
nated in Equation d35T > to the stars listed in Table 2. Below 
we show results of direct forward modeling; i.e., utilizing a 
known relationship for /* as a function of Rossby number. 
First, however, we wanted to investigate whether or not a 
single monotonic relationship for /*(Ro) could produce mass 
loss rates that were even remotely close to the measured val- 
ues. Thus, we produced trial grids of models in which /* was 
treated as a free parameter. For each star, we varied /* from 
10 to 1 and found the empirical value of the filling factor 
(femp) for which the modeled value of M matched the ob- 
served value given in Table 2. For the four binaries that have 
only systemic measurements of M we summed the model pre- 
dictions for each component and made a single comparison 
with the observations. 

Figure|9]shows the result of this process of "working back- 
wards" from the measured mass loss rates. The empirically 
constrained / emp values are plotte d against Rossby n umber, 
which is defined with (a) the simple Gunn et al.l (1 19981) ZAMS 
value for r c (Equation ( 1361 ) and (b) a gravity-modified ver- 
sion of t c that gives giants larger convective overturn times. 7 

7 We do not show / e mp for the test-case stars EV Lac or 89 Her, since 
no values in the range 10~ 5 -1 produced agreement with their observed mass 
loss rates. Extrapolating from the grid of modeled M values to the observed 
value would have required impossible values of / cm „ > 1000. We note, how- 
ever, that the F supergiant 89 Her "wants" to be in the unpopulated upper- 
right of the /<,(Ro) d iagram just like the F6 main sequence star HD 68456 
i Anderson et al. 2010). This may be relevant for deducing the relevant phys- 
ical processes in other F-type stars with T c ff es 6500 K. 



We varied the exponent in the gravity modification term and 
found that multiplying the ZAMS r e by (go/,g) ' 18 gives the 
narrowest distribution of /* versus Ro. The optimal exponent 
of 0.18 is very close to th e value of 0. 23 that w e found repro - 
duced the results of EH (1 19941) and iGondoinl d2005l 120071) . 
In Figure [9] we also show the same curves from Figure [7fb) 
that outline the measured range of filling factors. The lower 
envelope curve f m [ n , defined in Equation ( [381 , appears to be 
a good match to the gravity-modified empirical values / e mp- 
This provides circumstantial evidence that / m i„ is indeed an 
appropriate proxy for the filling factor of open flux tubes as a 
function of Rossby number. 

We now put aside the empirical estimates for the filling fac- 
tor and use only Equation (1381 for /* = f m [ n in the remainder 
of this paper. Table 3 shows some of the predicted properties 
of stellar coronae and winds for the 47 stars in our database. 
There were only seven stars for which Equation d38l gave a 
filling factor below the adopted "floor" value of /basal = 10~ 4 ; 
we replaced f m m with /b asa i in those cases. Table 3 also gives 
fk.TR, the coronal heat flux deposited at the TR for each star. 
It may be useful to use this to predict the X-ray flux associated 
with open-field regions on these stars, but we should note that 
the closed-field regions (which we do not model) are likely to 
dominate the observed X-ray emission. We also list the vari- 
ous components of Equation (f35j so that the contributions of 
gas pressure and wave pressure can be assessed (see below). 

Figure[10]compares the theoretical and measured mass loss 
rates with one another as a function of L» . For the four bi- 
nary systems listed in Table 2 with combined A+B mass loss 
rates, we separated the measured value into two pieces us- 
ing the modeled M ratio for the two components. (This was 
done for this figure only because the measured rates are shown 
as a function of a single star's luminosity.) For stars with 
only vsinz rotation period estimates, we used Equation (l40l to 
compute M for the central plotting symbol and the entries in 
Table 3, and we recomputed M for the lower and upper limits 
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Table 2. Cool Star Mass Loss Data 



Name 


T cff (K) 


log g 


M,/M e 


R./R® 


U/Lq 


P m (d) 


[Fe/H] 


-log(M/[M yr -1 ]) 


Ref. 


Sun 


5770 


4.44 


1 


1 


1 


25.3 





13.5-13.7 


— 


a Cen A (G2 V) 


5886 


4.31 


1.105 


1.224 


1.622 


29 


+0.197 


13.40 (A+B) 


33, 34, 23, 7 


a Cen B (KO V) 


5473 


4.54 


0.934 


0.863 


0.603 


36.2 


+0.230 


13.40 (A+B) 


33,34,35,7 


70 Oph A (KO V) 


5300 


4.52 


0.89 


0.86 


0.53 


19.7 


+0.040 


11.70 (A+B) 


14, 34, 23 


70 Oph B (K5 V) 


4390 


4.65 


0.73 


0.67 


0.15 


34 


+0.040 


11.70 (A+B) 


14, 34, 36 


e Eri (K4.5 V) 


5094 


4.60 


0.83 


0.754 


0.345 


1 1.7 


-0.097 


12.22 


18, 34, 23,7 


61CygA(K5V) 


4425 


4.63 


0.69 


0.665 


0.153 


35.4 


-0.193 


14.00 


22, 34, 23, 7 


e Ind (K5 V) 


4635 


4.54 


0.70 


0.745 


0.231 


22 


-0.088 


14.00 


37, 34, 23, 7 


36 Oph A (K5 V) 


5135 


4.54 


0.602 


0.69 


0.299 


20.3 


-0.206 


12.52 (A+B) 


19, 34, 7 


36 Oph B (K5 V) 


5103 


4.58 


0.486 


0.59 


0.213 


22.9 


-0.195 


12.52 (A+B) 


19, 34, 7 


i Boo A (G8 V) 


5551 


4.57 


0.86 


0.801 


0.550 


6.2 


-0.122 


13.00 (A+B) 


12, 34, 38, 7 


£ Boo B (K4 V) 


4350 


4.80 


0.70 


0.550 


0.0977 


1 1.5 


-0.122 


13.00 (A+B) 


12, 34, 38, 7 


61 Vir (G5 V) 


5560 


4.39 


0.946 


0.972 


0.804 


29 


-0.002 


14.22 


39, 34, 36, 7 


<5 Eri (KO IV) 


5025 


3.75 


1.122 


2.33 


3.185 


55.3 


+0.069 


13.10 


40, 34, 23, 7 


a Boo (Kl. 5 III) 


4290 


1.76 


1.10 


23 


170 


447 


-0.53 


9.60 


41,42,7 


a Tau (K5 III) 


3898 


1.33 


1.5 


44 


394 


648 


-0.180 


10.83 


43,44 


7 Dra (K5 III) 


3985 


1 .53 


3.0 


49 


535 


557 


-0.150 


1 1 .06 


43,44 


HR 6902 (G9 lib) 


4900 


1.99 


3.86 


33 


566 


220 


+0.430 


10.68 


45, 46, 47 


/3 And (MO III) 


3742 


1 .55 


1.5 


34 


204 


188 


-0.04 


10.19 


48, 49, 7 


UMi (K4 III) 


4040 


1.27 


1.3 


44 


475 


1030 


-0.26 


10.01 


48, 50,51,7 


fj. UMa (MO III) 


3700 


0.69 


1.5 


92 


1430 


488 


0.00 


8.92 


48, 49, 7 


a TrA (K2 II-IH) 


4150 


1 .50 


23.7 


143 


5500 


888 


-0.06 


9.77 


52,53,7 


A Vel (K4 Ib-II) 


3820 


0.64 


7.0 


210 


8510 


1250 


+0.23 


8.52 


54, 55, 7 


BD+01 3070 (RGB) 


5130 


2.70 


0.749 


6.4 


25.6 


50.9 


-1.85 


8.76 


56,57 


BD +05 3098 (RGB) 


4930 


2.00 


0.746 


14.3 


109 


109 


-2.40 


8.91 


56, 57 


BD +09 2574 (RGB) 


4860 


2.10 


0.753 


12.8 


82.5 


204 


-1.95 


8.94 


56.57 


BD +09 2870 (RGB) 


4600 


1.40 


0.864 


30.7 


381 


235 


-2.37 


8.43 


56,57 


BD +10 2495 (RGB) 


4920 


2.12 


0.723 


12.3 


80.0 


156 


-1.83 


9.06 


56, 57 


BD+12 2547 (AGB) 


4610 


1.50 


0.780 


26.0 


275 


265 


-0.72 


8.15 


56, 57 


BD+17 3248(RHB) 


5250 


2.21 


0.458 


8.80 


53.1 


64.8 


-2.02 


8.95 


56, 57 


BD +18 2757(AGB) 


4840 


1.43 


0.446 


21.3 


225 


127 


-2.19 


8.24 


56, 57 


BD+18 2976 (RGB) 


4550 


1 .30 


0.769 


32.5 


408 


248 


-2.40 


7.65 


56, 57 


BD-03 5215 (RHB) 


5420 


2.60 


0.884 


7.80 


47.4 


42.5 


-1.66 


8.68 


56, 57 


HD 083212 (RGB) 


4550 


1.40 


0.663 


26.9 


280 


146 


-1.49 


8.14 


56, 57 


HD 101063 (SGB) 


5070 


3.40 


1.19 


3.60 


7.73 


28.6 


-1.13 


9.56 


56, 57 


HD 107752 (AGB) 


4750 


1.70 


0.838 


21.4 


210 


185 


-2.88 


8.24 


56, 57 


HD 1 10885 (RHB) 


5330 


2.50 


0.757 


8.10 


47.8 


39.3 


-1.44 


8.66 


56, 57 


HD 111721 (RGB) 


5080 


2.35 


0.460 


7.50 


33.8 


74.5 


-1.26 


9.20 


56, 57 


HD 115444 (RGB) 


4750 


1.62 


0.584 


19.6 


176 


169 


-2.77 


8.56 


56, 57 


HD 119516 (RHB) 


5440 


2.37 


0.626 


8.60 


58.4 


39.8 


-2.50 


8.38 


56, 57 


HD 121135 (AGB) 


4925 


1.90 


0.789 


16.5 


144 


76.3 


-1.57 


8.23 


56, 57 


HD 122956 (RGB) 


4600 


1.35 


0.447 


23.4 


221 


130 


-1.78 


8.03 


56, 57 


HD 135148 (RGB) 


4275 


0.80 


0.239 


32.2 


312 


164 


-1.90 


7.83 


56, 57 


HD 195636 (RHB) 


5370 


2.17 


0.442 


9.10 


62.1 


16.9 


-2.83 


8.05 


56. 57 


EV Lac (M3.5 V) 


3168 


4.80 


0.315 


0.369 


0.0124 


4.38 


-0.200 


13.70 


26, 27, 34 23 


V Hya (N6, AGB) 


2160 


-0.89 


4.20 


945 


17540 


576 


+0.10 


5.12 


58, 59, 60 


89 Her (F2 lb) 


6550 


0.60 


0.61 


64.8 


6970 


143 


-0.41 


8.00 


61, 62, 63, 7 



(l)-(32) See Table 1, (33) Porto de Mello et al. (2008), (34) Wood et al. (2005a), (35) DeWarf et al. (2010), (36) Baliunas et al. (1996), (37) Janson et al. (2009), 
(38) Wood & Linsky (2010), (39) Vogt et al. (2010), (40) Hekker & Aerts (2010), (41) Schroder & Cuntz (2007), (42) Carney et al. (2008), (43) Robinson et al. (1998), 

(44) Cayrel de Strobel et al. (2001), (45) Kirsch et al. (2001), (46) Griffin (1988), (47) Marshall (1996), (48) Judge & Stencel (1991), (49) Massarotti et al. (2008), 
(50) Tarrant et al. (2008), (51) de Medeiros & Mayor (1999), (52) Ayres et al. (2007), (53) Harper et al. (1995), (54) Carpenter et al. (1999), (55) Setiawan et al. (2004), 
(56) Dupree et al. (2009), (57) Cortes et al. (2009), (58) Bergeat & Chevallier (2005), (59) Lambert et al. (1986), (60) Barnbaum et al. (1995), (61) Sargent & Osmer 

(1969), (62) Danziger & Faber (1972), (63) Stasinka et al. (2006). 



given by Equation ( P+TT i to obtain the error bars. Figure [TO^b) 
shows a comparison with th e semi- empirical scaling law pro- 
posed bv lSchroder & Cuntzl (|2005l) . with 

/ r e ff \ / , #o 



M sc = r/- 



1 + - 



(42) 



M* V4000KV V 4300g, 
where L», and M* are assumed to be in solar units. For 
this plot we computed the normalization constant r\ such that 
the average modeled mass loss rate would equal the average 
measured mass loss rate for all 47 stars. Averages were taken 
using the logarithm of M so that all stars would contribute to 
the average comparably. We found 77 = 8.5 x 10~ 14 yr" 1 , 
which is within the error bars of the lSchroder & Cuntzl (2005) 
value. 



Overall, our "standard model" (i.e., Equation (135b with 
cto = 0.5, h = 0.5, and 9 = 1/3) appears to match the mea- 
sured mass loss rates reasonably well. We emphasize that this 
model does not contain any arbitrary r\ normalization factors. 
For the three test-case stars at the bottom of Table 2, however, 
our model does not do as well. The model underpredicts the 
mass loss from the dM flare star EV Lac by at least four orders 
of magnitude, and it also fails for the F supergiant 89 Her by 
a slightly smaller amount. For EV Lac and other flare-active 
M dwarfs, it is possible that coronal mas s ejections and other 
episodic sources of energy (iMullanll 19961) could be responsi- 
ble for the bulk of the observed mass loss. For the carbon star 
V Hya, the reasonably good agreement between the model 
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TABLE 3 

Theoretical Wind Properties of Cool Stars 



Name 


Ro 


lOS /min 


B t (G) 


logoff TR 


Afhnt I 'MtvAA 

1T -*I1UI/ 1T -*LU1U 


log Ma tr 


-logfM/TMmyr'l) 


Sun 


1.960 


-2.996 


1513.05 


6.14 


19.73 


-2.42 


13.44 


a Cen A (G2 V) 


2.074 


-3.078 


1308.79 


6.17 


11.43 


-2.17 


13.21 


a Cen B (KO V) 


1.755 


-2.835 


1545.97 


5.62 


51.73 


-2.93 


14.09 


70 Oph A (KO V) 


0.996 


-2.025 


1666.45 


5.99 


194.1 


-3.25 


13.43 


70 Oph B (K5 V) 


1.027 


-2.067 


2130.57 


4.55 


417.3 


-4.50 


15.17 


e Eri (K4.5 V) 


0.519 


-1.174 


1832.80 


5.97 


977.4 


-3.94 


13.31 


61 Cyg A (K5 V) 


1.107 


-2.173 


2145.92 


4.59 


246.1 


-4.43 


15.15 


e Ind (K5 V) 


0.755 


-1.646 


1941.88 


5.15 


551.7 


-4.22 


14.25 


36 Oph A (K5 V) 


0.923 


-1.918 


1927.79 


5.67 


228.1 


-3.63 


13.84 


36 Oph B (K5 V) 


1.021 


-2.059 


1973.54 


5.51 


173.1 


-3.66 


14.16 


£ Boo A (G8 V) 


0.381 


-0.853 


1788.75 


6.69 


1335 


-3.59 


12.42 


| Boo B (K4 V) 


0.340 


-0.755 


2620.92 


4.83 


5964 


-5.37 


14.68 


61 Vir (G5 V) 


1.765 


-2.843 


1524.08 


5.99 


27.96 


-2.63 


13.56 


8 Eri (KO IV) 


1.844 


-2.907 


1077.96 


5.77 


12.04 


-2.18 


12.72 


a Boo (K1.5III) 


4.217 


-4.000 


411.54 


5.23 


0.57 


-0.31 


10.33 


a Tau (K5 III) 


4.183 


-4.000 


270.35 


5.17 


0.75 


0.11 


9.88 


7 Dra (K5 III) 


4.097 


-4.000 


301.38 


5.24 


0.85 


-0.12 


9.99 


HR 6902 (G9 lib) 


3.157 


-3.692 


287.82 


6.11 


1.28 


0.13 


9.86 


P And (MO III) 


1.229 


-2.321 


311.31 


5.61 


0.47 


-0.85 


8.86 


UMi (K4 III) 


6.960 


-4.000 


262.45 


5.32 


0.85 


0.28 


9.72 


H UMa (MO III) 


2.182 


-3.152 


165.83 


5.81 


0.63 


0.68 


7.93 


a TrA (K2 II-III) 


7.010 


-4.000 


270.87 


5.55 


1.25 


-0.13 


9.32 


A Vel (K4 Ib-H) 


5.808 


-4.000 


136.66 


5.67 


2.15 


1.05 


8.47 


BD +01 3070 (RGB) 


1.122 


-2.192 


792.36 


6.49 


2.36 


-1.48 


10.14 


BD +05 3098 (RGB) 


1.600 


-2.700 


553.26 


6.37 


0.50 


-0.63 


9.08 


BD +09 2574 (RGB) 


3.001 


-3.618 


632.02 


5.74 


0.49 


-0.62 


10.17 


BD +09 2870 (RGB) 


2.249 


-3.196 


476.21 


6.00 


0.37 


-0.17 


8.68 


BD +10 2495 (RGB) 


2.390 


-3.284 


602.67 


6.01 


0.51 


-0.62 


9.83 


BD +12 2547 (AGB) 


2.657 


-3.439 


378.07 


5.98 


0.60 


0.10 


9.20 


BD +17 3248 (RHB) 


1.259 


-2.355 


493.93 


6.81 


0.65 


-0.53 


9.05 


BD+18 2757 (AGB) 


1.400 


-2.508 


384.43 


6.55 


0.32 


-0.03 


8.05 


BD +18 2976 (RGB) 


2.218 


-3.175 


464.10 


5.96 


0.32 


-0.12 


8.54 


BD -03 5215 (RHB) 


1.095 


-2.157 


583.84 


7.06 


1.91 


-0.90 


9.36 


HD 083212 (RGB) 


1.361 


-2.467 


464.86 


6.20 


0.23 


-0.44 


8.07 


HD 101063 (SGB) 


0.813 


-1.744 


1312.86 


6.26 


40.43 


-2.65 


11.31 


HD 107752 (AGB) 


2.171 


-3.145 


524.60 


6.08 


0.42 


-0.37 


9.05 


HD 110885 (RHB) 


0.909 


-1.897 


575.05 


7.06 


2.02 


-0.97 


9.17 


HD 111721 (RGB) 


1.379 


-2.486 


609.98 


6.42 


0.61 


-0.91 


9.64 


HD 115444 (RGB) 


1.919 


-2.965 


491.45 


6.14 


0.34 


-0.31 


8.82 


HD 119516 (RHB) 


0.945 


-1.950 


483.08 


7.24 


1.29 


-0.62 


8.80 


HD 121135 (AGB) 


1.071 


-2.126 


502.44 


6.68 


0.60 


-0.68 


8.47 


HD 122956 (RGB) 


1.218 


-2.309 


444.03 


6.27 


0.19 


-0.38 


7.88 


HD 135148 (RGB) 


1.033 


-2.075 


394.38 


5.99 


0.07 


-0.25 


6.98 


HD 195636 (RHB) 


0.350 


-0.779 


435.75 


7.64 


3.38 


-0.88 


7.90 


EV Lac (M3.5 V) 


0.0706 


-0.313 


3005.25 


2.72 


2641 


-8.15 


17.78 


V Hya (N6, AGB) 


0.609 


-1.367 


12.39 


6.53 


1.32 


1.99 


4.34 


89 Her (F2 lb) 


27.12 


-4.000 


43.59 


1.92 


0.05 


-1.11 


11.15 



and measurements is probably a coincidence, since our model 
does not include the dusty radiative transfer or strong radial 
pulsations that are likely to be important for AGB stars. 

It is interesting to highlight the case of the moderately 
rotating K dwarfs e Eri, 70 Oph, and 36 Oph, which 
iHolzwarth & Jardinel d2007l) found to have anomalously high 
mass loss rates. They concluded that the observed magnetic 
fluxes for these stars were insufficient to produce their dense 
outflows. For these stars, our modeled mass loss rates tended 
to be about a factor of 10-20 below the measured values. 
However, these models were computed using / m ; n from Equa- 
tion J38T >. The measured values of /» given in Table 1 for these 
stars are larger than their corresponding / m ; n values by factors 
ranging from 3 to 20. If instead these values were used, our 
modeled mass loss rates woul d be in bette r agreem ent with 
the astrospheric observations of IWood et alj d2005al) . 

We also developed a statistical measure of how well a given 
model agrees with the measured database of mass loss rates. 



We defined a straightforward least-squares parameter 
1 N 

X 2 = ]jJ2 [logM,(model) - logM,(obs)] 2 (43) 
i=i 

where the total number of comparisons (N = 40) excludes the 
final three test cases in Table 2 and counts each of the four 
A+B binaries as one. When x 2 < 1, then (on average) the 
modeled and measured mass loss rates are within an order of 
magnitude of one another. Table 4 summarizes the results, in- 
cluding comparisons with other published empirical prescrip- 
tions. The rj normalization factors for each of these scaling 
laws were computed similarly as the factor in Equation d42l 
above. Note that our standard model appe ars to be a sig- 
nifica nt i mprovement over both the popular iReimersi (I1975L 
119771) and lSchroder & Cuntzl (120051) scalings. 

To further explore the proposed model, we varied some of 
the modeling parameters described in Section [3] Varying the 
TR filling factor exponent 6 did not have much of an effect 
on x 2 . However, varying the flux height scaling factor h did 
change \ 2 significantly. We found that a larger value of h w 
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FIG. 9. — Empirical / em p filling factors computed to match measured mass 
loss rates. Rossby numbers were computed in two ways: (a) directly from 
Equation (36), and (b) multiplying t c from Equation j36t by (g0/g) 018 . 
Symbol shading is proportional to log(L»/LQ), and the curves are the same 
as those in Figure [3b). The thick orange bar denotes the Sun's empirical 
range of values (computed from the small variation in M). The open orange 
circle denotes V Hya, the only one of the three test cases (at bottom of Table 
2) that gave a realistic solution for / era p. 



TABLE 4 
Mass Loss "Goodness of Fit" 



Model 


x 2 


This paper (standard model) 


0.650 


This paper (ZAMS t c ) 


1.575 


This paper (h = 0.25) 


0.794 


This paper (h = 1) 


0.564 


This paper (h = 3) 


0.504 


This paper (0 = 0.2) 


0.620 


This paper (0 = 0.5) 


0.707 


This paper (all [Fe/H] = 0) 


0.647 


This paper (M = M ho , +M co ld) 


0.703 


This paper (a from Trampedach & Stein 201 1) 


0.770 


Reimers (1975, 1977) 


1.260 


Mullan (1978), Equation (4a) 


3.768 


Nieuwenhuijzen & de Jager (1990) 


2.356 


Catelan (2000), Equation (Al) 


1.924 


Schroder & Cuntz (2005) 


1.131 
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FIG. 10. — Comparison of modeled (open circles) and measured (blue 
crosses) mass loss rates for the stars in Table 2, plotted as logM versus stellar 
luminosity. Panels show (a) the stan dard model developed in t his paper, and 
(b) the empirical scaling relation of Schroder & Cuntz (2005). The Sun is 
shown as a filled black circle, and the three test-case stars from the bottom of 
Table 2 are shown in orange. Vertical er ror bars in (a) correspond to models 
computed for the P rot range of Equation )4U . 

3 gives much better agreement with the measured mass loss 
rates than does the standard value of h = 0.5 (see Table 4). We 
decided not to adopt this larger value, though, because it falls 
well outside the range of empirically determined h values for 
the Sun's corona. 

We also tried removing some of the imposed complexity of 
the standard model to see if simpler assumptions would give 
adequate results. Removing the gravity-dependent modifica- 
tion factor of (go/g)° 18 from the definition of the convective 
overturn time resulted in significantly poorer agreement with 
the data (i.e., more than double the \ 2 of the standard model). 
We explored the importance of metallicity by replacing the 
published [Fe/H] by purely solar values ([Fe/H] = 0). This ac- 
tually improved the value of \ 2 from the standard model, but 
only by < 1%. Removing the exponential factor in Equation 
d35l l gave a slightly higher value of x 2 (8% larger than the 
standard model). 

We also n oted that the theoretical phot o spheric Alfven wave 
fluxes from Musielak & Ulmschneiderl d2002al) exhibited a 
strong dependence on the convective mixing length parameter 
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FIG. 1 1. — Illustrations of the relative importance of "hot" versus "cold" 
mass loss mechanisms, (a) Ratio of hot to cold modeled values for M 
(open circles) compared with the modified ratio M^ot ex P(~ 4^a trV^coM 
(red filled circles), (b) Mach number at the TR, Ma.tr = "tr/Va,tr- The Sun 
is shown as a thicker black circle in (a) and a filled black circle in (b). 

(i.e., Fa* oc a 21 ). Thus, instead of simply assuming a = 2 as 
in the standard model, we creat ed a linear regression fit to th e 
tabulated simulation results of Trampedach & Steinl (1201 lb . 
who found empirical values of a between 1 .6 and 2.2 depend- 
ing on r e ff , log g, and M» . We used the following approximate 
fit ' 

a T s ~ 1.91- + + — — (44) 



6181 K + 5.58 ' 19. 1M, 



and did not allow ctxs to be less than 1.6 or greater than 2.2. 
Thus, we multiplied the value of Fa* from Equation (|7]i by 
a factor of (axs/2) 2 ' 1 . The predicted mass loss rates for the 
Table 2 stars had about an 18% higher value of \ 2 than the 
standard model, so we did not pursue this mixing length pre- 
scription any further. 

For additional context ab out t he hot and cold mass loss 
models described in Section 13.31 Figure QT| shows the ratio 
Mhot/Mcoid for the 47 modeled stars as well as the TR Mach 
number M A TR = Mtr/Va,tr- It is clear that the dwarf stars 
are dominated by hot coronae, and the stellar wind outflow is 
still negligibly small at the coronal base. However, as the lu- 
minosity exceeds ~ 50 L Q for the giant stars, the hot coronal 
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FIG. 12. — Plasma number density at the TR plo tted as a function of stellar 
rotation period. Symbols are the same as in Figure [Tol 

contribution goes away and the acceleration becomes domi- 
nated by wave pressure. 

Figure Q~2] examines how the plasma number density at the 
transition region, wtr = Ptr/'«h, depends on stellar rotation. 
iHolzwarth & Jardind (120071) assumed «jr oc OP 6 oc P~® t 6 (see 
also llvanova & T aam 2003). Although our models do not fol- 
low a single universal relation for both giants and dwarfs, the 
proposed scaling (or one slightly steeper) may be appropri- 
ate for certain sub-populations of stars. For the dwarf stars 
with well-determined rotation periods, it is interesting that the 
Sun's computed value of «tr is larger than that of stars having 
higher magnetic activity. Equation ( 122b shows that the depen- 
dence on filling factor is weak (i.e., about f% 19 for 9=1 /3), 
so the variation comes mostly from the other stellar parame- 
ters. This seems to stand in contrast with other observational 
determinations of cor onal electron densities, where n tends to 
increase with activity (Giidel 2004). However, X-ray determi- 
nations of number density are probably dominated by closed- 
field active regions, which are not necessarily correlated with 
the regions driving the stellar wind. 

5.3. Predictions for Idealized Stellar Parameters 

In addition to the above comparisons with the individual 
stars of Table 2, we also created some purely theoretical sets 
of stellar models and computed M for t hem. We began with 
the ZAMS model parameters given by iGirardi et al.l d2000), 
and we assumed solar metallicity for a range of constant rota- 
tion rates. This gave rise to a two-dimensional grid of models 
(varying T e & and P rot ) for main sequence stars. Because the 
modeled stars are all high-gravity dwarfs, we used only Equa- 
tion (l36l > for r c , in combination with Equation (1381 for / m j n as 
a function of Rossby number. 

Figure Q~3] shows the resulting mass loss rates as a func- 
tion of T e ff and P 10t . If we had not utilized a basal "floor" on 
/*, we would have predicted a steep drop-off in m ass loss 
for r eff > 7000 K, at which the iGunnetalJ (119981) expres- 
sion for r e decreases rapidly. However, because of the floor, 
there appear to be reasonably strong mass loss rates up to 
the point at which subsurface convection zones disappear at 
?eff ^ 9000 K. There is a slightly discontinuous dip in the 
predicted basal mass flux around T e ff w 6100 K that arises be- 
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FIG. 13. — Theoretical predictions of M for main sequence stars, plotted 
as a function of T c ff. Differently colored curves show a range of assumed 
rotation periods (see labels for values). The Sun is indicated by an open 
circle. 

cause of the iteration for 1Z, pjR, and £?tr- If the calculation 
of these quantities is halted after only one iteration, the final 
value of M varies more smoothly as a function of r e ff. We plan 
to utilize a more self-consistent non-WKB model of Alfven 
wave reflection in future versions of this work. 

The mass loss rates shown in Figure[l3]are a lmos t all due to 
the hot coronal processes discussed in Section [3~Tl Thus, it is 
possible to simplify the components of Equation d26l i in order 
to obtain an approximate scaling relation for M that is rea- 
sonable for these main sequence stellar models. Ignoring the 
weakest dependences on some stellar parameters (i.e., factors 
with exponents less than or equal to 1 /7), we found 



M 



10- 10 M Q /yr 



R 



16/7 



•Pa* 



12/7 



-2/7 



10 9 erg cm 2 s 



(4+30)/7 



(45) 



which reproduces the curves in Figure [13] to within about an 
order of magnitude. Despite the fact that this scaling formula 
is relatively easy to apply, we do not recommend its use in 
stellar evolution or population synthesis calculations. Once 
stars leave the main sequence, Equation d45l > is no longer a 
good approximation. 

We also computed a time-dependent mass loss rate for the 
evolutionary track of a star having M* = 1M Q . There is 
evidence that the wind from the "young Sun" was signifi- 
cantly denser than it is today, and this more energetic out- 
flow may have been important to early planetary evolution 
(e.g..lWoodl2006 : Gudel 20071 ISterenborg et al.l201 lllSuzukil 
1201 ll). We used the BaSTI ev olutionary track plotted in Fig- 
ure[T](Pietrinferni et al. 2004) for the time variation of R* and 
L*. We grafted on a model of rotational evolution for a solar- 
mass star from Figure 6(a) of Denissen kov et al.l (120 1 Oh . For 
late ages (f > 100 Myr, or logf > 8), this model has approx- 
imately P rc 



oc t 



0.54 



Such an age scaling is well within the 

') ob- 
tained from young solar analogs (e.g., iBarnesI 120031; iGiidell 
2007). 
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FIG. 14. — Theoretical predictions of stellar wind properties for an evolv- 
ing solar-mass star, (a) Base-10 logarithms of luminosity (green dot-dashed 
curve), Rossby number (black solid curve), and / m j n (blue dashed curve) plot- 
ted as a function of age in years, (b) Our standard mo del for M (black solid 
curve), compared with the Schroder & Cuntz 12005) scaling (blue dashed 
curve) and an ideal power-law t 1 decline with increasing age (green dotted 
curve). A model and observations of CTTS are shown for comparison (red 
error bars and dot-dashed curve). 

FigurefPfta) shows how the luminosity and two dimension- 
less parameters related to the rotational dynamo (Ro and / m i n ) 
vary as a function of age for this model. Note that prior to 
about t ps 70 Myr the Rossby number is small enough that 
the filling factor appears to be saturated near its maximum 
assumed value of 0.5. At very late times, when the star be- 
gins to ascend the red giant branch, the Rossby number de- 
creases again because of the increase in t c with decreasing 
r e ff and gravity. We utilized the (g©/g) ' 18 correction factor 
when computing r c , but it was relatively unimportant until the 
star left the main sequence. 

Figure fT4T b) gives our prediction for the age variation of 
the a solar-type star's mass loss rate. For ages between about 
t fa 0.2 and 7 Gyr the decrease in mass loss appears to be 
fit approximately by a power law, with M oc t" 1 . This is a 
significantly shallower age depe ndence than the f" 2 decline 
suggested by IWood et alJ (120021) on the basis of astrosphere 
measurements. We note that if the rotation period was the 
only variable to change with time, Equation d45l l would give 
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something like M oc P~^ A (for /* = / m ; n and 9=1 /3). Thus, 
a more rapid increase of f ro t with age — such a s the t° 85 de- 
pendence in the solar-mass rotational model of lLandin et alj 
(2010) — would give rise to a steepe r age-M relationship more 
similar to that of iWood et all (120021) . 

For comparison, Fig ure \l4\ b) also shows that the 
ISchroder & Cuntj j2()05) scaling law predicts a much smaller 
range of mass loss variation for the young Sun than does 
the p resent model. We al so show a mod el (Cranmerl l2008l 
2009) and measurements (Hartiga n et al.|[l995l) for classical 
T Tauri stars (CTTS) at the youngest ages. It is clear that 
for t < 10 Myr some additional physical processes must be 
included (e.g., accretion-driven turbulence on the stellar sur- 
face) to successfully predict mass loss rates. 

6. DISCUSSION AND CONCLUSIONS 

The primary aim of this paper was to develop a new genera- 
tion of physically motivated models of the winds of cool main 
sequence stars and evolved giants. These models follow the 
production of MHD turbulent motions from subsurface con- 
vection zones to their eventual dissipation and escape through 
the stellar wind. The magnetic activity of these stars is taken 
into account by extending standard age-activity-rotation indi- 
cators to include the evolution of the filling factor of strong 
magnetic fields in stellar photospheres. The winds of G and 
K dwarf stars tend to be driven by gas pressure from hot 
coronae, whereas the cooler outflows of red giants are sup- 
ported mainly by Alfven wave pressure. We tested our model 
of combined "hot" and "cold" winds by comparing with the 
observed mass loss rates of 47 stars, and we found that this 
model produces better agreement with the data than do pub- 
lished scaling laws. We also made predictions for the para- 
metric dependence of M on T e ff and rotation period for main 
sequence stars, and on age for a one solar mass evolutionary 
track. 

The eventual goal of this project is to provide a straightfor- 
ward algorithm for predicting the mass loss rates of cool stars 
for use in calculations of stellar evolution and population syn- 
thesis. A brief stand-alone subroutine called BOREAS has 
been developed to implement the model described in this pa- 
per. This code is written in the Interactive Data Language 
(IDL) 8 and it is included with this paper as online-only mate- 
rial. This code is also provided, with updates as needed, on 
the first author's web page. 9 Packaged with the code itself are 
data files that allow the user to reproduce many of the results 
shown in Section[5] 

In order to further test the conjecture that Alfven waves 
and turbulence drive cool-star winds, the models need to 
be expanded from the simple scaling laws of Section [3] 
to fully self-consistent solutions of the mass, momentum, 
and energy conservation equations along open flux tubes. 
Modeling the full radial dependence of density, temper- 
ature, magnetic field strength, and outflow speed would 
eliminate our reliance on approximat e factors like h and 
9. We believe that th e models of ICranmer et al.1 (120071) 
for the solar wind, and ICranmerl (|2008|) for T Tauri stars, 
can be extended straightforwardly and applied to other 
types of stars. However, there are many other approaches 
to producing self-consistent and/or three-dimensional mod- 

8 IDL is published by ITT Visual Information Solutions. There are also 
several free implementations with compatible syntax, including the GNU 
Data Language (GDL) and the Perl Data Language (PDL). 

9 http://www.cfa. harvard.edu/~scranmer/ 



els that should be explored (e.g ., lAirapetian et al.1 
20 lOt iHolzwarth & Jardinel 12005b ISchriiver & Titli 



2000 



2005 



Falceta-Goncalves et al. | 120061; ISuzukil 120071; IVidotto et al.l 
2009; Cohen et alj|2009[ ICohedl2oTlir 

There are additional ways t hat o ur simplified models of 
coronal ener gy b alance (Section l3~TT l and wave-pressure driv- 
ing (Section [3!2l may be improved: 

1. Our standard assumption for the outflo w speed in a 
coronal wind was m„o = V esc . However. Hudgel (1 19921) 
found that many stars have significantly smaller termi- 
nal speeds. It should be poss i ble to use something like 
the ISchwadron & McCom as (2003) solar wind scaling 
law to estimate the peak tem perature in the corona, and 
thus apply the Parkej] (119581) theory of gas pressure ac- 
celeration to compute the wind speed. 



2. We assumed in Se ction 13.11 that rjR ~ R*. However, 
ISchroder & Cuntzl (120051) estimated that some low- 
gravity stars should exhibit "puffed up" chromospheres 
with a fractional extent given by the final term in paren- 
theses in Equation d42t . We applied this correction fac- 
tor to the modeled values of ttr and Mhot for the stars in 
Table 2. Doing so yielded significant differences from 
the standard model only for stars having M co id 3> Mhot, 
i.e., the combined model value of M was relatively un- 
changed in those cases. However, in general there may 
be other stars for which this kind of correction factor 
needs to be considered in more detail. 

3. We also assumed that the flux height scaling factor h 
took on a single constant value fo r all stars. It may 
be use ful to explore extending the ISchroder & Cunta 
(2005) idea of gravity-dependent spatial expansion to 
this parameter as well. The \ 2 results shown in Table 4 
suggest that a larger value of h sa 3 could be appropriate 
for many of the low-gravity stars in our observational 
database, whereas the range h « 0.5-1 is probably best 
for main sequence stars like the Sun. 

4. The flux of energy Fa* in kink/Alfven waves in the pho- 
tosphere may depe nd on other parameters that we have 
not considered. Mu sielak & Ulmschneiderl (l2002al) 
found that the flux is rather sensitive to B*/B eq in the 
photosphere, so departures from our assumed value of 
1.13 may give rise to significantly differe nt predictions. 
Also, Musielak & Ulmschneider (2002b) examined the 
sensitivity to metallicity and found that lower Z/Z Q 
tends to give lower values of Fa* for T e s < 6000 K. Pre- 
liminary tests showed that this effect does not strongly 
affect the mass loss rates derived in this paper. How- 
ever, a varying metallicity should also change other 
properties of the convection — including the effective 
mixing length a parameter — thus possibly making this 
effect more important. 

5. Instead of assuming a simple monotonic dependence 
of the open-field filling factor on Rossby number, it 
may be possible to construct realistic surface distri- 
butions of active regions for a given activity level 
and rotation period, and model the opening up of 
flux tubes by both stellar win d s and centrifugal forces 
( Mull an & Steinolfsonl ll983UMestel & Spruitl 119871: 
Jardinel 120041: IHolzwart h & Jardinel 120051 : ICohen et alj 
20091) . 
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To continue testing and refining these models, it is also im- 
portant to utilize the newest and most accurate measure- 
ments of stellar mass loss rates (see, e.g., [Schroder & Cuntzl 



2007; Willson 2009; Catelan 2009; Mauron & Josselin 201 1 



Vievtes et al. 2 0TT|) and magn etic fields dDonati & Landstreel 
2009tlVlemmings et alJl201 ll) . 



Finally, we emphasize that a complete description of late- 
type stellar winds requires the incorporation of other phys- 
ical processes besides Alfven waves and turbulence. The 
outer atmospheres of cool stars are also likel y to be pow- 
ered by acoustic or longitudinal MHD waves (ICuntzlll990t 
iBuchh olz et al. 1998 1), episodic flares o r coronal mass ejec- 
tions (lMullanirT 996; Aarn io et alJl2009l). and large-amplitud e 
pulsations dBowenlll988l;lde Jager etaTll 1 997t IWillsonl(200dl) . 
It is well known that radiative driving should not be neglected 
for AGB stars and red supergiants, an d it may be important 
for Cepheids (iNeilson & L ester 2008) and horizontal branch 



stars dVink & Ca ssisi 2002) as well. 
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APPENDIX 

NOTES ON STELLAR MAGNETIC FIELD MEASUREMENTS 

In this work we fo cus on observations of unpolarized spectral lines sensitive to Zeeman splitting by stellar magnetic fields 
(e.g., [Robinson 1980). The resulting "Zeeman broadened" line profiles are valuable probes of both the intensity-weighted mean 
absolute value of the field strength (i.e., B* = (7b|B|)/(/b), where 7b is the continuum intensity in the magnetic regions) and 
the intensity-weighted fraction of the visible stellar hemisphere that is covered by these fields (i.e., the filling factor /»). These 
detections are thus weighted towards the brightest regions of the stellar surface; i.e., plage or network regions. However, this 
technique allows detection of more topologically complex fields, and thus more comprehensive values of /* and B„, than does the 
use of circular polarization. The latter exhibits significant signal cancellation when there are multiple oppositely directed patches 
of magnetic field in the same resolution element. In many cases, however, only the disk-averaged magnetic flux density (B*j*) 
can be determined re liably from Zeeman broadened spectra and not the separate values of B* and /* (see also lRiiedi et al.lll~997l : 
I Anderson et~aT1l2010l) . 

Many deta ils of the observations of the stars discussed in Section|4]were given bv lSaar & Linskvl (fl^fl986h andHaai] (fl99(i[ 
1199 ll 1996a b, 2001). The approximate quality factors listed in Table 1 span the range from low (q = 1) to high (q = 4) relative 
confidence in the derived magnetic parameters. The values for q were assigned based on a combination of the following properties 
of the spectroscopic data and its magnetic analysis: 

1. Detections using lines with longer wavelengths and higher Lande g e ff factors are given higher q values. The ratio of the 
strength of the Zeeman effect to the nonmagnetic Doppler width is AAb / AAd, where AAb oc g e ff\ 2 is the Zeeman splitting 
amplitude and AAd oc A is the nonmagnetic Doppler width. Thus, the relative detectability of the effect increases as g e ff A. 

2. Spectra with higher signal-to-noise (S/N) ratios and higher spectral resolution (TZ = A/ A A) tend to have higher quality 
(see, e.g., Saar 1988), although a longer wavelength measurement can trump better 1Z. For example, the great er magnetic 
sensitivity at large A can lead to partial resolution of the individual Zeeman components dSaar & Linsky||1985h . 

3. The simultaneous analysis of larger numbers of lines (especially with higher g e ff) contributes to a good quality measurement 
dRiiedi et alJll997l) . It is additiona lly helpful f or the lines to be free of blends and for any rotational broadening to be small 
(i.e., vsku < 10 km s" 1 ); see also Saar ( 1988). 

4. Finally, the method and quality of line modeling — including the level of detail in the magnetic radiative transfer, the atmo- 
spheric models used, details of integration over the stellar disk, and the treatment of line blends — vary widely. This field 
has seen a continual improvement in the sophisticatio n of these models, with an attendant incr ease in our understanding of 
the measurement and systematic uncertainties (e.g.. lAnderson et aT1l2010t IShulyak et al.ll2.Q10h . 

We chose to exclude several p ublished observations from the list of stars given in Table 1 . We did not use the M dwarf data 
presented by iReiners et al.1 ([2009) because these stars all tended to sit in the saturated region of Figure[7](i.e., Ro/Ro Q < 0.1), 
and thus do not contribute to improvi ng our knowledge of th e rotation dependence of /*. We did not include the measured field of 
the F6 main sequence star HD 68456 (Ander son et al.l20 10) because we do not attempt to model the outflows of stars significantly 
hotter than the Sun. Also, its combined strong field and high Rossby number point to a possible transitio n to a different type 
of ma gnetic activity from that described by the standard cool-star age-activity-rotation relationship (see also Bo hm-Vitense et all 
2002)? 
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